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SHORT-TERM CLIMATIC FLUCTUATIONS 
FORCED BY THERMAL ANOMALIES 

The aim of this research is to study the response of the atmosphere 
to thermal anomalies using a low order spectral model. Thermal anomaly 
patterns may exist either in sea and land surface temperatures or in the 
tropospheric diabatic heating. 

A two-level, global, spectral model using pressure as a vertical 
coordinate has been developed. The system of equations describing the 
model is nonlinear and quasi-geostrophic (linear balance). Static 
stability is variable in the model. A moisture budget is calculated in 
the lower layer only. Convective adjustment is used to avoid super- 
critical temperature lapse rates. The mechanical forcing of topography 
is introduced as a vertical velocity at the lower boundary. Solar 
forcing is specified assuming a daily mean zenith angle. The differen- 
tial diabatic heating between land and sea is parameterized. On land- 
and sea-ice surfaces, a steady state thermal energy equation is solved 
to calculate the sui face temperature. On the oceans, the sea surface 
temperature is specified as the climatological average for January. The 
model is used to simulate the January, February and March circulations. 

Experiments are designed to study the response of the atmosphere to 
thermal anomalies at the lower boundary or in the midtroposphere. The 
**memory*’ in the atmosphere of such anomalies, after they have decayed, 
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is also studied. Three patterns of sea surface temperature anomalies 
are tested. The first pattern represents a cold anomaly in the North 
Pacific, the second a warm anomaly in the equatorial Pacific and the 
third pattern contains both of the two anomaly patterns acting together. 
The results suggest that the coupled pattern is the only one that pro- 
duces the type of geopotential anomalies associated with tie negative 
phase of the Southern Oscillation. In contrast to the results of linear 
models, warm sea surface temperature anomalies in the equatorial Pacific 
cannot produce such geopotential response on their own. In che case of 
this tropical anoiai'ly pattern the variance of temperature resulting from 
transient eddies tends to increase, whereas in the case of the coupled 
anomaly pattern the variance of temperature resulting from stationary 
eddies increases. This behavior suggests that with both anomalies 
acting together the atmosphere is inclined to produce quasi-permanent 
responses, such as blocking, in contrast with the other case, in which 
the transient activity increases. 

The mid-tropospheric anomaly is introduced as an easterly pro- 
pagating wave over the equatorial Pacific and over the Gulf of Bengal. 
The amplitude and memory of the response is larger than for the sea 
surface temperature case. The mid- tropospheric thermal anomalies show 
continuous large areas of long memory in the subtropical and middle 
latitude regions of the Northern Hemisphere. 
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1 . INTRODUCTION 


Climate is a subject of considerable importance to human lives. 
The impact of different climatic regimes on the economy, health and 
wealth of individuals and nations cannot be denied. Of the same impor- 
tance is the variability of climate. This variability has been ob- 
served, directly or indirectly, over a wide band of time scales ranging 
from seasons to thousands and millions years. Theories have been 
proposed in order to explain the phenomena associated with climatic 
variability on different time scales. These theories could be categor- 
ized into two main groups, depending on terrestrial or extraterrestrial 
forcing mechanisms. The terrestrial-forcing climate change theories 
have been concerned with continental drifts, mountain building, volcan- 
ism, variations of sea level and sea temperature, and variations in the 
carbon dioxide content of the atmosphere. The second category of 
theories has been concerned with the variations in the intensity of 
solar radiation and variations in the earth's orbit around the sun 
(Milankovitch, 1941). No completely acceptable explanation of climatic 
change has ever been presented. The problem will remain until we gain a 
complete understanding of the interactions between thf earth, atmo- 
sphere, oceans, and the solar and planetary system. Lorenz (1968) 
explored the possibility that the atmosphere may be an almost "intran- 
sitive system". In such a system the statistics taken over long, but 
finite, intervals differ considerably from one interval to another. If 
this is the case, there is a possibility that long-term climatic changes 
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may result from the nearly intransitive nature of the atmosphere, rather 
than from environmental changes. However, this theory applies to cases 
where environmental influences are not affected, in turn, by the system. 
The atmosphere by itself is not a system of this sort since the oceans, 
and also glacial ice areas, are affected by the atmospheric conditions. 

In this thesis we will be focusing on short-term climatic variabil- 
ity, namely the variability through a month or season. The response of 
the atmosphere to sea surface temperature anomalies (SSTA) and to tropo- 
spheric thermal anomalies will be explored. 

1.1 Observational Studies 

The sensitivity of the atmosphere to long-lasting perturbations in 
the lower boundary has been the subject of many investigations. SSTA 
and areas of ice coverage are the most important of these perturbations. 

Observational studies show possible te? econnections between SSTA 
and the atmospheric circulation. From analyses of three winter cases 
(1957-1958, 1963-1964 and 1965-1966), Bjerknes (1969) showed that posi- 
tive SSTAs, in the equatorial central and eastern Pacific Ocean in late 
fall were accompanied by an anomalous strength of the midlatitude west- 
erlies over the northeast Pacific. According to his hypothesis, these 
positive SSTAs increase the local vertical motions, thus enhancing the 
Hadley circulation, which in turn transports more than the average 
angular momentum to middle latitudes. Observations have shown that such 
wind anomalies are much more pronounced in the Northern Hemisphere than 
in the Southern Hemisphere. Nonlinear interactions of the velocity and 
thermal fields are more significant in the Northern Hemisphere midlati- 
tudes than in the Southern Hemisphere. The former has higher and 
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broader mountains than the latter, as well as larger differential diaba- 
tic heating between land and oceans. Canton Island data showed a maxi- 
mum rainfall. Thus the atmosphere has shown both local and remote re- 
sponse to equatorial SSTA. 

Namias (1969, 1980) studied the effects of midlatitude SSTAs in the 
central Pacific on the weather in the United States of America during 
the following season It was found that positive SSTAs in the central 
Pacific that occurred from 1961 through the winter of 1967-1968 enhanced 
and displaced the region of cyclogensis southward. This pattern caused 
a blocking ridge over the western United States. The ridge, with its 
maximum amplitude at 160* W, was accompanied by cold advection over the 
eastern United States, with mild temperatures prevailing over the cen- 
tral and western states. Cold SSTAs during the summer of 1978 in the 
central northern Pacific were followed by below-normal temperatures over 
all the United States in winter. The same studies showed that the 
response of atmospheric pressure to cold SSTAs would be a higher than 
normal pressure at the surface and aloft during the following fall. 
Uniform temperature anomalies over the entire United States usually are 
rare since the dominant long waves usually have smaller dimensions than 
the lateral extent of the United States. Thus the western and eastern 
parts of the country often display opposite temperature anomalies. 

Reiter (1978) confirmed the feedback between the Intertropical 
Convergence Zones (ITCZ), the trade winds and sea surface temperatures. 
That study further indicated that this feedback is only a part of a 
global atroo.^phere-ocean interaction, among many other factors. Analysis 
of two severe winters (1976-1977 and 1977-1978) in the eastern United 
States, by Harold (1980) indicated that transient eddy activity during 
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abnormal winters usually is not significantly different from such acti** 
vity during normal winters. The height and temperature anomalies were 
in phase, resulting in a structure equivalent to that of barotropic 
waves. Deviation of the jet stream from its normal position explained 
the observed precipitation patterns during the two winters. 

The sea surface temperatures over the North Atlantic and North 
Pacific have roughly the same variance as surface air tesq>erature on 
monthly, seasonal and annual time scales (Cayan, 1980). Wan-cheng ^ 
al . (1981) showed that increased (decreased) SSTA in the equatorial 
Pacific is associated with a delay by two months and possibly also with 
simultaneous, strengthening (weakening) of the upper-tropospheric west- 
erlies in the eastern Pacific subtropics. They did not find any associ- 
ation between SSTA and anomalies of the northward momentum transport in 
the subtropical upper tropospheric flow over the eastern Pacific. 

So far, the response of the atmosphere to SSTAs has been shown to 
be both local and remote. For the atmosphere to respond in such a 
manner, there must be planetary zonal or meridional oscillations that 
control and organize the phase and amplitude of such a response. Cor- 
relations between geopotential heights on a given pressure surface at 
widely separated points on the earth confirmed the presence of at least 
four patterns which might be indicative of standing oscillations in the 
planetary waves during the Nothern Hemisphere winter (Wallace and 
Gutzler, 1981). Examples are found over the North Atlantic, the North 
Pacific, polar temperate latitudes, and between the Pacific and North 
America. The North Atlantic oscillation was first identified by Walker 
and Bliss (1932) and confirmed by amoy studies using different indices 
for the oscillation (Van Loon and Rogers, 1978; Kutzbach, 1970; Wallace 


and uUtzler, 1^81). Positive values of the index are indicative of a 
strong Icelandic low, high pressure along 40®N, strong westerlies across 
the North Atlantic, below normal temperatures in the Greenland Labrador 
area as well as in the Middle East, and above-normal temperatures in the 
eastern United States and northwestern Europe. The North Pacific oscil- 
lation is an analogue of the North Atlantic oscillation (Walker and 
Bliss, 1932; Rogers 1981). Both the North Atlantic and North Pacific 
oscillations are characterized by north’south seesaws or standing oscil- 
lations in the sea level pressure with a node located near 50® latitude. 
The zoually symmetric seesaw (Lorenz, 1951) is basically the negative 
correlation between sea level pressure at polar and teipperate latitudes. 
This correlation also exists in the abser.ce of the North Atlantic and 
North Pacific oscillations. The last oscillation of this kind is the 
Pacific/North American pattern (Wallace and Gutzler, 1981), which essen- 
tially describes the blocking phenomena. Blocking high pressure is 
identified by the persistence of a quasi-stationary positive anomaly in 
the 500 mb geopotential height field. This anomaly persists for a 
period of at least 10 days (i.e., has a time scale longer ban the 
typically associ.ated with synoptic-scale variability) over at le**st one 
location (Dole, 1978). This phenomenon often happens during winter 
along the west coast of North America. The blocking is accompanied by a 
strong Aleutian low over the central Pacific and below normal geopoten- 
tial heights over the southeastern United States. 

To complete the discussion of the teleronnection patterns around 
the globe we must include two important oscillations that dominate the 
tropical and subtropical regions. The correlation in the surface pres- 
sure between eastern and western hemispheres at low latitudes was 
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recognized a long time ago in what is called the Southern Oscillation 
(Walker, 1923, 1924, 1928). It is defined as the tendency of the sea 
level pressure at stations in the Pacific (San Francisco, Tokyo and 
South America), anC of rainfall in India to increase, while pressure in 
the region of the Indian Ocean decreases. The positive index of the 
oscillation is recognized as negative sea ^evel pressure anomalies in 
the Indonesian belt with respect to positive sea level pressure anoma* 
lies in the Pacific belt. The existence of the Southern Oscillation was 
proved after Walker's discovery (Troup, 1565; Berlage, 1966; Horel and 
Wallace 1981). The Southern Oscillation has a period betv’een 3 and 6 
years, and it is not confined to tropical and subtropical latitudes, but 
it is linked to midlatitude systems and long waves in the westeriiis 
(Trenberth, 1976). The Southern Oscillation is a manifestation of a 
nearly global variation in circulation, temperature, clouds and pre- 
cipitation. Such a low frequency oscillation is linked to the fluctua- 
tions in sea level pressure in the equatorial Pacific. The ocean tem- 
perature is a triggering influence in developing the oscillation. When 
Walker discovered the Southern Oscillation in the twenties he did uot 
realize the source of "memory" for such a very low frequency phenomenon, 
since all the atmospheric fluctuations are of much shorter periods. The 
question remained without answer until the work by Bjerknes (1969), 
mentioned above, where he emphasiy.ed that the oscillation draws its 
memory from the long-lasting ocean temperature anomalies in the equa- 
torial Pacific. El Nino is a surface oceanic phenomenon recognized as 
warm water accumulated near the coast of Peru just before or after 
Christmas, '^his warm water interrupts cue normal upwolling. El Nino is 
preceded by strong southeast trades at the center of the Pacific ocean. 
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which intensify the subtropical gyre of the South Pacific, strengthen 
the South Equatorial Current, and increase the east-west slope of sea 
level by building up water in the western equatorial Pacific. As soon 
as the wind stress in the central Pacific relaxes, the water flows 
eastward. An El Nino event lasts for about 16 months (that is, two 
summers and one winter) (Wyrtki, 1975). 

The Walker circulation is a counterpart of the Southern Oscillation 
in the wind field. It is characterized by large scale rising motion 
over the warm zone in the west Pacific, near Indonesia, and sinking 
motion over regions of cold water in the east Pacific. The main forcing 
mechanism for the Walker circulation is the release of latent heat of 
condensation over the equatorial west Pacific. The release of latent 
heat enhances the vertical motion and results in upper-level divergence. 
The upper-level westerlies and the low-level easterlies over the Pacific 
ocean force a large scale subsidence over the cooler water of the south- 
east Pacific. The tropical circulation is not zonally synmietric due to 
the existence of such a circulation. Gill (1980), using a linear primi- 
tive equation model, speculated on the idea that Kelvin waves emanating 
from the Indonesian region would create an extensive region of easter- 
lies the east (which extend over tn.'^ whole equatorial Pacific), while 
planetary waves would create a more limited region of westerlies to the 
west (over the Indian Ocean). 

Many studies have shown possible teleconnections between the South- 
ern Oscillation, Walker Circulation, El Nino and tropospheric pressure 
patterns at different areas. Morel and Wallace (1981) identified the 
atmospheric pattern associated with both the local and remote response 
of the atmosphere during seven warm episodes, within the period 
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1951-1978, in the Pacific Ocean surface temperatures. T J mb geopo- 
tential was taken as an indication of the tropospheric thickness in the 
tropical region. Fluctuations in mean tropospheric temperatures and 200 
mb heights were shown to vary simultaneously with equatorial Pacific sea 
surface temperature fluctuations, not only in the Pacific sector, but at 
stations throughout the tropics. Negative SSTA in the northwestern and 
central Pacific at latitude 50°N occurred coincident with the east 
Pacific warm equatorial SSTA. The mid-tropospheric pattern associated 

with this sea surface temperature distribution showed below normal 

« 

heights in the North Pacific and southeastern United States and above 
normal heights over western Canada. The extratropical response was 
dependent on the southward shift of the raidlatitude jet stream to reach 
the lower latitudes. Fritz (1982) analyzed the same data of Morel and 
Wallace, but using the 700 mb average height difference between years of 
warm SSTA and years of cold anomalies. The analysis was done using a 
monthly average instead of a winter average. The pattern was most 
pronounced during January and February, rather than December. The 
stronger correlation during those months between the 700 mb height and 
sea surface temperature over the Pacific may be attributed to snow and 
ice distribution, which affect the atmospheric circulation. Analyses by 
Egger et (1981) of satellite data collected between February 1965 
and July 1973 revealed a minimum amount of cloudiness near the date line 
(180®) when the Southern Oscillation index was positive. At the same 
time, the wind maps showed an anomalous easterly flow. 

A positive Southern Oscillation index and a strong Walker circula- 
tion are associated with a southwest shift of the SoutJi Pacific 
Convergence Zone (SPCZ) from its position between 125 ®E, 0®S and 
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150®E> 23®S (Gruber, 1972). The ITCZ is displaced away from the equator 
into the Northern Hemisphere (Trenberth, 1976). Satellite cloud pic- 
tures show two tropical convergence zones in the west and central 
Pacific, south and north of the equator, respectively. The Indonesian 
continent and eastern Australian regions have a strong surge of pre- 
cipitation, while the area east of 160®E is dry. On the other hand, 
with a negative index of the Southern Oscillation, the SPCZ is pushed 
from its mean tropical position towards the equator and the ITCZ moves 
into a near equatorial position (Ramage, 1975). Cold water upwelling at 
the equator is replaced by positive SSTA in the equatorial central and 
eastern Pacific (El Nino). The equatorial dry zone east of 160®E ex- 
periences one of its prolonged precipitation surges. On satellite 
photographs one extended convergence zone marked by a cloud band near 
the equator extends across almost the entire width of the Pacific ocean. 

These tropical phenomena seem to have a teleconnection with mid- 
latitude and extratropical geopotential height fields. Reiter (1981) 
showed that low height anomalies dominate the 500 mb central North 
Pacific pattern for months, with precipitation maxima in the Line 
Islands, The 500 mb planetary wave anomalies in midlatitude can serve 
as triggers for tropical precipitation patterns, especially during 
January. Deep troughs over the central North Pacific during that month 
apparently have a tendency to push the ITCZ equatorward. A similar 
conclusion was reached by Van Loon and Rogers (1981). The negative 
phase of the Southern Oscillation is associated with a high zonal-mean 
poleward flux of sensible heat, and the 700 mb temperature and heights 
tend to be lower between 30°N and 60®N with a minimum value at 45°N. 
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The interannual variability of precipitation in the tropics has 
been attributed to the variability of the ITCZ (Kraus, 1977; Hastenrath, 
1976; Hastenrath and Heller, 1977). A recent study by Stoeokenius 
(1981) showed that this hypothesis failed to explain many inportant 
rainfall anomalies in the tropics. In other words the zonally averaged 
Hadley circulation does not seem to be able to account for most of the 
observed interannual rainfall variations, and hence of general circula- 
tion variations, in the tropics. East-west circulations like the 
Southern Oscillation and the Walker Circulation have an important in- 
fluence . 

1.2 Theoretical Investigations 

The quasi-stationary ultra-long waves (zonal wave number 1-4) are 
of major importance in determining the large scale characteristics of 
the atmosphere. The interannual variability of their phase and ampli- 
tude is an indication of the change of the behavior of the atmosphere 
and, therefore, of regional climate. In the Northern Hemisphere wave 
number 2 is the dominant energy peak, while in the Southern Hemisphere 
it is wave number 1 (Holton, 1975). The development of atmospheric 
ultra-long waves is due to lower boundary forcing, namely the earth's 
surface. Differential heating between land and sea together with the 
distribution of topography are the two main mechanisms that force the 
ultra-long waves (Sankar-Rao, 1965; Sankar-Rao and Saltzman, 1969). 
This also was confirmed by Ashe (1979) using a spectral model which 
solved a time average asymmetric flow induced by topography and diabatic 
heating. When the forcing of topography was considered alone, it was 
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found that most of the features were displaced from their true posi- 
tions. In middle latitudes the latent heating (calculated from pre- 
cipitation averages) tends to be in phase with the sensible heating. 
The opposite is true in the subtropics, particularly in desert regions, 
where strong positive departures of sensible heating are balanced by 
strong negative departures of latent heating. Tung and Lindzen (1979) 
stressed the importance of the non-linear interactions, and forcing by 
the topography and diabatic heating in simulating the major circulation 
features. During a normal year the two kinds of forcing act out of 
phase. This is because the major north-south mountains in the Northern 
Hemisphere (e.g., the Rockies and Alps) lie along the eastern side of 
the Pacific and Atlantic oceans, respectively. The major warm ocean 
currents, like the Kuroshio and Gulf Stream, flow in the western side of 
the Pacific and Atlantic. With this normal forcing structure, the two 
mechanisms seem to act against each other. The zonal wind speed at 500 
mb is larger than the phase speed of the ultra-long waves and the pat- 
tern propagates eastward with the zonal flow. Since the topographic 
distributions are quasi-permanent, if the diabatic heating due to land 
and sea distribution was altered, then we would expect the two kinds of 
forcing to act in harmony to reinforce, instead of cancelling, each 
other. If the wind circulation in the lower atmosphere is such that the 
phase speed of the wave is reduced to zero, a nearly resonant state may 
be reached by the amplifying planetary waves, causing a blocking pheno- 
menon. In all the previous studies the mean zonal flow was specified; 
not calculated by the models. 

We notice in the preceding discussion that the fundamental feature 
of the forcing of stationary waves in the atmosphere is that the forcing 
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is largely implicit, depending on the wind circulation in the atmosphere 
itself. Charney and DeVore (1979) and Hart (1979) stressed this fact. 
In their theory of "multiplo equilibrium states", blocking was found to 
be one of many equilibrium states that can be reached by the system 
itself. Roads (1980), using a nonlinear quasi geostrophic model in 
which orography and zonal heating are calculated, has found multiple 
equilibrium states. The states are either near a resonance point in the 
system or at the zonal fot''ing equilibrium point. The multiplicity of 
the solution occurs because the heat transport and mountain torque are 
nonlinear functions of the zonal stream function. More analysis was 
done by Charney et al. (1981) using a barotropic nonlinear channel model 
incorporating observed zonal topography. Multiple stationary equilibria 
were obtained. From the observed 34 blocking events of the 500 mb geo- 
potential, only 19 could be explained as one or another of the equilib- 
rium states, whereas 5 events could be explained due to forcing or 
geometry. What is not explained is the localized position of the ridges 
and troughs and mechanisms of transition to and from the blocking con- 
figuration. In other words, the theory could not prove whether or not 
there is (are) a certain mechanism(s) responsible for forcing, or 
ending, any or all of these multiple solutions, including blocking. 

Linear models with steady state solutions have been used to iden- 
tify certain responses of the atmosphere to large scale forcing. There 
is some agreement between the linear model results and the observed 
anomalies. However, the linear theory is questionable for such types of 
studies due to the strong nonlinearity of the atmosphere itself. It is 
doubtful that the positions of the actual extrema caused by large scale 
anomalies could be accurately predicted. Egger (1977) tested the effect 
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of tropical SSTA in both the eastern Pacific and the Atlantic using a 
linear two-level hemispheric model. In the Pacific case the response of 
the midlatitude atmosphere to the tropical anomalies was negligible. 
This result was not supported by observations. However, the tropical 
Atlantic anomalies showed a larger response by the midlatitude atmo- 
sphere. 

The sphericity of the earth could be an important factor in deter- 
mining the correct response of the atmosphere to the large scale anoma- 
lies, in particular the propagation of Rossby waves on a sphere. The 
results suggest that the nature of the response is quite different 
between the different zonal sectors of the globe. Webster (1981) pre- 
sented the argument that, in low latitudes, where the convective acti- 
vity and the diabatic heating are the dominant factors in the thermal 
energy equation, the response of the atmosphere to the anomalies seems 
to be of a direct nature. In other words, the warm SSTA in the tropical 
oceans enhances the vertical velocity, which in turn increases the 
poleward momentum flux in the upper layer of thi troposphere and 
strengthens the nidlatitude westerlies. In this case the teleconnec- 
tions between low and midlatitudes are strong and clear. On the other 
hand advection is of primary importance in the thermo-dynamic energy 
equation in middle latitudes, and the anomalies cause an indirect cir- 
culation (rising of cold air and sinking of warm air) rather than a 
direct circulation as in the low latitude case. Therefore, the local 
response to the midlatitude anomalies is smaller than that for the low 
latitude case. The teleconnection between the midlatitude atmosphere 
and midlatitude SSTAs depends upon the nature of the atmospheric circu- 
lation. An example is whether the anomalies are in a region of rising 
or sinking motion. 


Subtropical SSTA has an effect on the atmosphere in middle and high 
latitudes since the response of the atmosphere to the tropical heating 
can have a substantial influence on the middle and high latitudes, 
provided that part of the heating is in the belt of westerly winds 
(Opsteegh and Van Den Dool, 1980). In this case, modes with meridional 
components to the group velocity vector will be excited which will allow 
a substantial response at higher latitudes (Webster, 1982). The propa' 
gation of the long wave lengths as a result of subtropical thermal 
heating was shown to be polewards as well as eastwards, while shorter 
wavelengths were trapped equatorward of the poleward flank of the jet 
(Hoskins and Karoly, 1981). 

The question of the generally larger response of the atmosphere to 
the tropical anomalies than to midlatitude anomalies still has no satis- 
factory answer. The previous linear models (Opsteegh and Van Den Dool, 
1980; Webster, 1981, 1982; Hoskins and Karoly, 1981) are constrained by 
the assumption of zonal synanetry in the basic state. Observational 
studies of general circulation statistics indicate that quasi-stationary 
modes in the extratropical regions and the tropics provide significant 
longitudinal perturbations in the mean state of the atmosphere. Such 
quasi-stationary perturbations probably will alter the character of both 
the teleconnection propagation and the local form of the atmospheric 
response . 

The importance of nonlinear aspects of the atmosphere has been 
demonstrated by many investigators. The amplitude of baroclinic waves, 
and ^o a lesser extent the amplitude of their transports of heat and 
momentum, are sensitive to the barotropic component in the basic flow 
(Simmons and Hoskins, 1980). Some of this sensitivity has been traced 
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to distinctly barotropic effects, which arise during the nonlinear 
evolution of the waves. Frederiksen and Sawford (1981) used a spherical 
equivalent barotropic m- ^ to study the linear and non-linear solutions 
for the effects of topography in forcing the stationary eddy field. 
There are considerable differences between the two types of solutions, 
which arise largely from the fact that linear solutions are resonant or 
near resonant, whereas the nonlinear solutions are not. Comparison with 
observed fields show nonlinear effects to be most important at low 
latitudes. 

1.3 General Circulation Models 

Apart from the simplified linear theory discussed before, many 
efforts have been devoted towards finding the sensitivity of the atmo- 
sphere to environmental anomalies using large scale general circulation 
models (GCMs). Several important studies of this type were the ex- 
periments done by Rowntree (1972, 1976b). He used a nine-level hemi- 
spheric model of the atmosphere developed at the Geophysical Fluid 
Dynamics Laboratory (GFDL) to test Bjerknes hypothesis that fluctuations 
of ocean temperatures in the tropical east Pacific are responsible for 
major variations in the position and intensity of the Aleutian surface 
low. The experiments were done with both warm and cold SSTA with a 
maximum difference of 3.5“C. The results of the experiments showed that 
tropical ocean temperatures have important effects on the model atmo- 
sphere, in the tropical and extratropical regions. The warm ocean 
induced a surface low over the area of maximum temperature, and the 
associated low level convergence and ascent gave increased tropical 
rainfall over the central and east Pacific, with decreases over the west 
Pacific and parts of South America. An upper-tropospheric subtropical 
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jet stream was generated with persistent trough-development to the north 
over the midlatitude east Pacific. The Aleutian low shifted east and 
deepened, confirming Bjerknes hypothesis. In such a hemispheric model 
with boundary conditions imposing symmetry about the equator, it was 
argued that the equatorial wall in Rowntree's model was responsible for 
the anomalous precipitation calculated in the tropical region (Ramage 
and Murakami, 1973; Ramage, 1975). However, the experiments using a 
global 5 layer GCM (United Kingdom meteorological office) (Corby et al. , 
1977) gave results consistent with the hemispheric experiment. 
Experiments with SSTA over the Atlantic agree with those done over the 
Pacific. Decreases of surface pressure were obtained over the warm 
tropical water and extended north to at least bO^N over the eastern 
Atlantic (Rowntree, 1976b). In many experiments it was noticed that 
warm anomalies gave more significant atmospheric response than cold 
anomalies (Houghton ^ al . , 1974). 

Many experiments were done for testing the effect of middle lati- 
tude and extratrupical SSTA. In most of these experiments it was dif- 
ficult to detect the signal above the noise level, but compositing of 
several experiments with different anomalies suggests a tendency for 
pressure falls on the equatorward side of the anomalies (Rowntree, 
1979). In other experiments testing North Pacific SSTA, the results 
suggest that a statisfically significant signal can be detected in the 
immediate vicinity of large ocean surface temperature anomalies, but 
that teleconnections of this signal downstream over the United States 
can be identified only with a small statistical confidence (Chervin 
et al . , 1976). As a matter of fact, extremely cold or warm extra- 


tropical SSTAs, were needed in both the Pacific and Atlantic Oceans to 
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detect a signal with a high confidence level (Kutzbach et aJL • , 1977; 
Chervin eX al . , 1980). In the last series of experiments by Chervin et 
al. (1980), experiments were done assuming a dipole-type anomaly distri- 
bution with warm anomalies in the west of the North Pacific and cold 
SSTA to the east. This pattern may not resemble the observed pattern 
(Horel and Wallace, 1980). At the same time the atmospheric response 
may be reduced since the cold (warm) anomalies are under a region of 
rising (sinking) air. 

Simpson and Downey (1975) used the NCAR model to find the sen- 
sitivity of the atmosphere to SSTA in the middle latitudes in the 
Southern Hemisphere. They used two patterns of SSTA in addition to that 
used in the control run. The SSTAs were associated with completely 
different atmospheric circulations. The atmosphere responded more to 
one pattern than the other. This may suggest that the response of the 
atmosphere to SSTA may depend upon the initial atmospheric circulation 
associated with the location of the anomalies. While many numerical 
models failed to support Namias’ (1969, 1980) hypothesis, the results of 
a recent experiment using the GLAS (Goddard Laboratory for Atmospheric 
Science) model (Shukla and Bangaru, 19811 supported his hypothesis. The 
SSTA and atmospheric circulation for this experiment were similar to 
those observed during January 1977. It was shown that a cold SSTA over 
the Pacific produces a strong southward flow over the United States and 
colder temperatures in eastern Canada and the United States, as was 
observed during winter 1977. Contrary to the previous numerical exj^er- 
iments, the results indicated that the SSTA over the Pacific can pro<luce 
a significant downstream effect over the continental United States. The 
different results between the NCAR experiment (Kutzbach ^ > 1977; 
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Chervin ^ al. , 1976, 1980) and the GLAS experiment may be attributed to 
the spatial structure of the SSTA. 

The more significant response of the atmospheric models to tropical 
SSTAs than to midlatitude and extratropical SSTAs still has no com- 
pletely natisfactory explanation. It is believed that both the physical 
and dynamical characteristics of the interaction between the ocean and 
atmosphere should be explored further. Kraus (1972) claimed that, via 
an enhanced air-sea interaction, the role of SSTA should be of most 
significance in the tropics for the following reasons: 

1) Perturbations induced by changes in the zonal wind stress will 
have a much faster response in the tropics due to the higher group 
velocity of oceanic Rossby waves at lower latitudes. 

2) Local upwelling should produce larger surface temperature 
variations because of the shallower thermocline in the warm sea 
areas. 

3) Surface temperature variations in the tropics can produce 
larger variations in the heat supply than is possible in higher 
(colder) latitudes. 

If ocean-atmosphere feedbacks are not considered, the middle lati- 
tude SSTA could be inefficient in significantly changing the global 
diabatic heating. The atmosphere over these areas behaves in a highly 
nonlinear fashion and is affected by transient eddies and baroclinic 
instability. Ashe (1978) discussed an indirect offset of middle lati- 
tude SSTA. The anomalies would change the static stability parameter 
near the surface, which in turn would affect the baroclinic instability 
of the atmosphere. The advection of unstable eddies of this kind could 
cause a significant effect on the longwave pattern, and this large scale 
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change could feed back to affect the diabetic heating itself (Ashe, 
1978). Cyclogensis at middle latitudes is created by land and sea 
surface temperature contrast, the zonal wind, the mountain chains and 
the latent heat released. So the relation between middle latitude SSTA 
and the atmosphere is more complicated than a simple direct response of 
the atmosphere to the anomalies. 

The complex and nonlinear interactions between different variables 
describing the ocean-atmosphere circulations are an important factor in 
obtaining significant results from different experiments. A change in 
the lower boundary surface temperature may affect the cloud distribu- 
tion, which again has a feedback to the atmospheric circulation. The 
results of a simple, coupled air-sea model (Lau, 1979) indicated that 
sea surface temperature distribution, which is strongly controlled by 
the oceanic upwelling, is the primary factor in determining the location 
and movement of the tropical rainbelt. The strongest convective ac- 
tivity in the ITCZ, however, depends mainly on the moisture supply from 
horizontal convergence and the static stability of the lower atmosphere, 
but do not necessarily coincide with the occurrence of maximum sea 
surface temperatures. The increase of sea surface temperature, while 
increasing the evaporation, condensation and precipitation, was found to 
be associated with a decrease in the global cloudiness and relative 
humidity (Schneider and Dickinson, 1974; Roads, 1978). Sellers (1976) 
suggested that a possible reason for these results was that an increase 
in temperature may correspond to an increase in convective cloudiness, a 
decrease in stratus cloudiness, and a decrease in total cloud area, 
along with an increase in condensation from increased evaporation. 
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Stephens and Webster (1981) concluded that surface teaq;>erature and ver> 
tical temperature structure is acutely sensitive to the cloud aaraunt and 
cloud form. Low clouds cool the surface and higl clouds warm it. A 
study by Wecherald and Manabe (1980) showed that if the solar constant 
was increased, clouds diminish in the upper and middle troposphere at 
most latitudes and increase in the lower stratosphere and near the 
earth's surface. The area mean change in the net incoming solar radia* 
tion is coriipensated by the outgoing terrestrial radiation at the top of 
the atmosphere. Owing to this compensation, the change of cloud cover 
hai> a relatively minor effect on the sensitivity of the area anan cli- 
mate. Meanwhile another study by Shukla and Sud (1981), rsing the GLAS 
model to investigate the effect of cloud radiation feedback, showed a 
significant change in the simulated large scale circulations of the 
model. 

Ramage (1977) suggested that interpretation of the precipitation 
patterns from numerical models must be made wxth caution because most 
numerical models are non- interactive; the sea surface temperature is 
prescribed and does not change ae a result of the atmospheric conditions 
above it. Even with interactive models, the problem still exists with 
convective pai...^eterization. The most complex parasieterizatlon schemes 
still cannot resemble the cosqplicated nature of convection, especially 
in the tropics. The results may give a realistic picture of the actual 
atmosphere. As an exasiple of this problem, Ramage tested an esipirical 
formula given by Rowtree (1976a). The empirical formula correlates 
precipitation with SSTA in the Pacific region east of the date line and 
from the equator to 20^N. In this region the area-averaged sea surface 
temperature increased 1.7°C between December 1971 and December 1972 . 
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According to the empirical formula of Rowntree's model, the rainfall 
should have increased 58 mm/day. But the observations showed that 
rainfall decreased by 3.4 mm/day. 

1.4 Scope of the Present Study 

The present research has two main phases. The first is the devel- 
opment of a global spectral model. The second is to examine theoretical 
arguments using such a model. We intend to gain more understanding of 
the dynamical response of the atmosphere to possible existing patterns 
of S^ A and mid-tropospheric thermal anomalies. The points of investi- 
gation are: 

1) The atmospheric response to SSTA at different locations. For 
e iample, how does the atmosphere respond to a cold SSTA in the 
North Pacific, a warm SSTA in the equatorial eastern Pacific, or to 
both at the same time? The response of the atmosphere to negative 
SSTA in the North Pacific and positive SSTA in the equatorial 
region was not examined before. A recent study by Horel and 
Wallace (1981) suggested that this pattern of the SSTA has a signi- 
ficant influence on Northern Hemisphere weather. 

7 } The effect of mid-tropospheric thermal anomalies which may 
t ;ist as a result of anomalous release of latent heat. Atmospheric 
whermal anomalies are more likely to occur as propagating waves 
which resemble the atmospheric motion in a certain region. We will 
be concerned with thermal anomalies in the tropical region which 
are triggered by the interannual variability of the location of the 
ITCZ. 

3) The atmosphere may have a different memory of the thermal 
anomalies at different locations. In other words, if the anomalies 
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were removed there would be different recognition of such an effect 
by the atmosphere, depending on their location and the associated 
atmospheric circulation. 

The experiments are designed to examine both the local and remote 
atmospheric response to SSTA or mid-tropospheric thermal anomalies. The 
analysis includes the horizontal wind, geopotential height, temperature, 
vertical motion, and evaporation fields. The average tropospheric 
amplitudes of the transports and variances resulting from stationary and 
transient eddy components also are analysed. 

For this purpose a low order spectral global model in linear bal- 
ance was developed. Nonlinear interactions, topographic forcing and 
sea-land differential heating were included. There are no assumptions 
of zonal symmetry or specification of the zonal flow. This type of 
medium complexity model, in which the dynamical and physical processes 
are simplified witfout contradicting important observations of the 
atmosphere, provide a further insight into the results of the compli- 
cated general circulation models, which usually have more degrees of 
freedom. The model developed has an advantage over steady state linear 
models or non-steady state linear models in that the time dependent 
nonlinear interaction of belting and motion can be simulated. On the 
other hand, the simplicity of the model relative to large scale general 
circulation models gives the advantage of an easier interpretation of 
the model’s results. 

The economic computer runs of medium complexity models have the 
advantage of examining more hypotheses and running for longer periods of 
integration without excessive cost. In the present study the model was 
used for a 3 month simulation of the winter circulation for the control 
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and anomaly cases. The completion of 30 days integration takes about 50 
seconds using the NCAR CRAY I machine. The model could be used eco- 
nomically for extended periods of integration. Kraus and Lorenz (1966) 
used a model by Lorenz (1960) for very long periods of integration to 
study the variability and predictability of the atmosphere. Models of 
such design have been used for climate sensitivity studies. Salmon and 
Hendershott (1976) used a quasi-geostrophic model to study the effect of 
SSTA on the atmospheric circulation. Kikuchi (1979) used a model 
similar to the one presented here^ except that in the present model we 
assume a variable static stability and a hydrologic cycle in the lower 
layer for studying the influence of orography and land-sea distribution 
on the Northern Hemisphere winter circulation. Held and Suerz (1978) 
developed a two level primitive equation model which was used to examine 
the response of the atmosphere to perturbations of the solar constant 
(Held et al. , 1981). 

Chapter 2 presents the model description and formulation. In 
Chapter 3 the simulated January climate of the model is compared with 
observations. Experiment design and results are explained in Chapters 4 


and 5. 


2. P TWO- LEVEL GENERAL CIRCULATION MODEL 


2. 1 Introduction 

A two level global, spectral model using pressure as a vertical 
coordinate has been designed. The system of equations retains the 
nonlinear interactions between the dependent variables. The system is 
quasi-geostrophic in linear balance form (Lorenz, 1960). The mechani- 
cal effects of graphy have been introduced in the form of a lower 
boundary vertical v locity. The differential diabatic heating due to 
the distribution of land and sea also are parameterized. Both orography 
and differential heating between land and sea are important mechanisms 
for producing a correct amplitude and phase cl the middle latitude 
ultra-long waves (wave numbers 1-4) (Sankar-Rau and Saltzman, 1969; 
Bates, 1977). The solar radiation budget is specified as a function of 
latitude and time. Longwave emissivities also are included. Static 
stability is variable in the model. This condition is necessary for 
conserving the sum of energies of the model. Also, it is important for 
our purpose of sensitivity studies. The model uses a moisture budget 
equation in its lower layer only, with a convective adjustment to pre- 
vent super-critical lapse rates. Under reversible adiabatic processes, 
the equations conserve the sum of kinetic energy and available potential 
energy. 

In the horizontal domain a rhomboidal spectral truncation is as- 
sumed (truncated at zonal wave number 9) (Fig. 2.1). The dependent 
variables are expanded in terms of their harmonic coefficients with 
spherical harmonics as basis functions. 
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2.2 Basic Equations 

The follov'ing system of equations forms the basic structure of the 
model. The velocity field at each level is split into its rotational 
and divergent parts. For a certain level, £, the system of equations 
representing the model's atmosphere is given by 


Vjg = kxVilijj + Vx^, (2.1) 

the vorticity equation 
a 

* 3T " * <V,- 


the thermodynamic energy equation 


fPof 

^ = -J(4.,,e,)-vx,.ve,-u,^ 8T Mi;] c " (^.3) 


the thermal wind equation 
'■o'* 

9p 

the Continuity equation 
9u)/i 

— ^ + V^Y = 0 

dp ^ H 

and the moisture budget equation (at the lower level only) 
|a=-V(Vjq)*E-P^MS^)j. 


(2.4) 


(2.5) 


( 2 . 6 ) 


In (2.1) to (2.6), V = (u,v) is the horizontal velocity, u) the 
”<»rtical pressure velocity, tj< the stream function, X the velocity poten- 
tial, 6 the potential temperature, and q the water vapor mixing ratio. 
(Fj^) and (F^) are the horizontal and vertical friction respectively. 
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f is the Coriolis parameter; K = where R is the gas constant for 
dry air, and is the specific heat at constant pressure, p^ is the 
pressure at the lower boundary (1000 mb). (W^^) and (W^) are the hori- 
zontal and vertical temperature diffusion, respectively. is included 
in tie diabatic heating rate, Q, after parameterizing the vertical heat 

fluxes. (S. )- is the horizontal diffusion of moisture. Q is the sum of 
h 1 

diabatic heating rates per unit mass. E is the evaporation rate from 
the surface to the lower layer and is the precipitation rate. 

The diffusion and diabatic heating functions need to be parameter- 
ized in terms of the dependent variables X£> and q in 
order to form, with the necessary boundary conditions, a closed system 
of equations. 

The two levels represented in the model are 750 mb (£=1) and 250 mb 
(£=3). The average values of the modeled variables are calculated at 
the lower and upper layers. An intermediate level, 500 mb (£=2) also is 
used for calculation of the vertically averaged conditions of the 
model *s atmosphere. The lower boundary is the 1000 mb level (i^=0) (Fig. 


2 . 2 ). 

The upper layer is heated by short- and long-wave radiation, by the 
**ateral diffusion of heat, and by the heat released by convective ad- 
justment. The lower layer is heated by short and longwave radiation, 
lateral diffusion of heat, seisible heat flux from the surface, and 
latent heat release due to condensation, and is cooled by the heat 
transferred upward by the convective adjustment. Evaporation from the 
surface is the only source of moisture to the lower layer. Condensed 
water is assumed to fail as precipitation, without any re-evaporation 


into the lower layer. 
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a.® 

1 

o 

n mb 

^3* ^3* ^3 

250 mb 

“2 

500 mb 

>i'i. Xj^» q 

750 mb 

0) 


S 

1000 mb 


Schematic representatioD of the vertical structure of 
the model. 
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2,3 Boundary Conditions 

In the horizontal resolution all wave number modes up to wave 
number 9 have been retained, in order to maintain the interaction be- 
tween tbe two hemispheres. A wall at the equator was not assumed. 

The boundary condition for the vertical velocity in a pressure- 
coordinate system, u), at the top of the model's atmosphere is 
U)4 = 0. 


The mechanical influence of orography is introduced in the model as a 


vertical velocity at the lower boundary (1000 mb), lo ^0, Rather, 

s 

U) = V *Vp , 
s ~s 

where is the rotational pait of the 7S0 mb wind velocity. The con- 
tribution from the divergent part is at least one order of magnitude 
less, hence is neglected, p^ is the surface pressure due to orography. 
Applying the continuity equation (2.5) at different levels, the vertical 
pressure velocities at the different levels are given by 

(2.7) 

^ ( 2 . 8 ) 

V2(2x3+Xj), (2.9) 

and u)^ = -Ap V2(x3+Xj), (2.10) 

where Ap is the difference between the upper and lower levels = 500 mb. 
We define x^» the velocity potential at the surface, such that 

( 2 . 11 ) 


u>2 = -ApV^x^, 




u». 


= -ApV^x^, 


Then 

= 2(v^-vn) 


( 2 . 12 ) 
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whi-re n is the normalized pressure = Pg/Pgi Pq “ 2Ap 
from (2.10) and (2.11), and we obtain 

X, = X, * Xj (213) 

For computing the surface pressure, p , the continental elevations, 

S 

smoothed over 5“ latitude x 5® longitude "squares" are used (Berkofsky 

and Bertoni, 1955). The conversion from elevation heights to pressure 

.3 

is done assuming an atmosphere with a constant lapse rate of 6.5 x 10 
®K m‘^ 

The low order truncation used in the model is considered as a 
further filter to satisfy the quasi-geostrophic approximation, where the 
vercical velocity should be three orders of magnitude less than the 
horizontal wind (Haltiner, 1971). Fig. 2.3 shows the smoothed orography 
distribution used in the model. 

2.4 Horizontal Diffusion 

The horizontal diffusion terms are not necessary for coaiputational 
stability when using a spectral model. They are required to inhibit 
spurious growth of amplitude at scales close to the point )f truncation 
due to spectral blocking (Puri and Bourke, 1974). At a level i the 
dix.fusion terms are given by 




(Fh)l = kh(V^(V^*l) - 2-^ ). 

(2.14) 


(2.15) 

and (Sjj)g = kj^ V^q. 

(2.16) 


2 *1 

where kj^ = 1.0 x 10 ^ m sec (Phillips, 1956) is the lateral eddy 
diffusion coefficient, and a is the radius of the earth. The last term 
on the right side of (2.14) is due to the effect of a spherical earth. 




32 


2.5 Vertical Diffusion 




2.5.1 Parameterization of Frictional Dissipation 

The two assumptions used for parameterizing the frictional dis~ 
sipation are as follows (Lorenz, 1961): 

a) Surface frictional drag is proportional to the flow in the 
surface layer; 

b) friction between the two layers is proportional to the 
differen'-i! between the flow of the two layers. 


(F„) = - g 


3p 


( 2 . 17 ) 


where r^ is tLe rotation stress due to vertical diffusion at 
level £. 




(2.18) 


^''2 ■ 


(2.19) 


The rotational stress at the top of the model's atmosphere, r^, is 
assumed to be equal to zero. The coefficients of friction at the under- 
lying surface and the surface separating the layers are denoted by 
and -kj, respectively. 

r. = + ^ k (2.20) 

0 g s ^s 

and ''d (2.21) 

where is the surface stream function extrapolated from the values at 
levels 250 mb and 750 mb, assuming (j« is linear with respect to height. 
‘>>ibstituting (2.20) and (2.21) into (2.18) and (2.19), we obtain 
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and 


(F,)^ = -k, (4.3 - ♦,), 


(F,) = kj Vk(»3 - *,) - k^ Fk*, 
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( 2 . 22 ) 


(2 23) 


where = pPj^ (^)^ and I'g = MPg > P small scale vertical 

eddy stress coefficient, and and p^ are the air densities at level 2 
and at the surface respectively. The two coefficients (Kikuchi, 1969 
and Charney, 1959) are given by 

k = 4 X 10 ^ sec \ 
s 

and = 5 X 10 ^ sec 

2.5.2 Surface Sens i ble and Latent Heat Fluxes 

The heat and moisture fluxes are calculated in the same way as was 

described by Bourke et ajL. (1977). The sensible heat flux at the lower 

layer is calculated using the bulk formula 

Q = p c IV I (T - T ), (2.24) 

^s ^s p d ‘ si g ^ ^ 

where T is the surface air temperature at anemometer level, T is 
a a 

extrapolated fron* the temperature values at 250 mb and 750 mb to level 
1000 mb, c^ is the drag coefficient, and T^ is the surface temperature, 
which is determined from the surface thermal energy balance. Over 
water, the sea surface temperatures are specified. 

The evaporation from the surface to the lower layer is given by 

Tk* ' amount of latent heat lost from the surface due lo evaporation is 


Qe = LE, 


,2.26) 


where h is a measure ol the relative huimdilv in th^ laver. The wet- 
s 

ness parameter, GW, is equal to unity on water surfaces* Over land 
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surfaces 


PAO£lS 


GW = 0.25. 

The drag coeffrcent is assumed to be constant over land and oceans. 
Over land c^ = 0.004 and for the ocean case c^ = 0.001 (Bourke ^ al. , 
1977). is the saturation mixing ratio at 1000 mb. 


2.6 Radiation Budget 


2.6.1 Solar Radiation 

The incoming solar radiation at the top of the model's atmosphere 
is calculated as a function of daily mean zenith angle (Wetherald and 
Manabe, 1972). Diurnal variation of the solar energy is excluded. 
The mean zenith angle z is given by 

cosz = sin$ sinfi + (cos<|> ccs6 sin H^)/H^, (2.27) 

where 0 is the latitude angle, 6 is the declination angle, and is the 
hour angle given by 

= cos ^ (-taru|> tan6) ^2 28) 

6 = 23.45 sin27T (|^) • 

N is the number of days measured from day 0 at OOZ 1 Janu_<.y. 

The incoming solar radiation at the top of the atmosphere is given 
by 

S = S H /n, (2.29) 

00 

where 


S = 



0-6 
0 - 6 


< n 
2 



(2.30) 
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2 

is the solar constant, assumed as 1400 W/m . Recent measurements of 

solar irradiance from earth orbiting satellites (Smith ^ al. , 1982) 

2 

give an average value about 1375 W/ro . 

a and a are the instantaneous and mean distances of the earth from 
s 0 

the sun, respectively, 

^ = 1 + .01676 sin2n (|^) • (2.31) 

o 

The solar radiation absorbed by the atmosphere (Saltzman, 1967), S^, is 
given by 

= x(l-r^)S, (2.32) 

where x is the absorptivity of the atmosphere, assumed to be constant = 

0.26 in this study. 

r^ is the albedo of the atmosphere. This parameter was calculated 
using a formula by Vernekar (1975), which includes the mean zonal amount 
of clouds. 

r = (a + pc)c (2.33) 

a 

p is a constant equal to 0.38, a is a function of latitude, c is the 

amount of low and medium clouds in tenths of sky cover. Although the 

model ha^ no explicit modulation of the clouds, they are implicitly 

included through the atmospheric albedo and the solar energy budget. 

The net solar energy absorbed by the earth's surface is given by 

= (1-x) (1-r ) (1-r ) S, 
s as 

where r^ is the January zonal-average albedo of the earth's surface 
(oceans are not included). The surface albedos are categorized as areas 
of permanent ice (albedo ~ 0.8), partial snow or ice cover (Albedo =0.4 
to 0.6), deserts or patchy snow in middle and low latitudes (Albedo = 
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0,2 to 0.3), dense forests (Albedo = 0.10 to 0.15). The values of 
different paraineters used for the January solar radiation calculations 
are given in Table (2.1). 

The above formulae give a global average planetary albedo = 34%. 
Stephens ^ (1981), using satellite observations, estimated the 

global average planetary albedo for January tc be 31%. Fig. 2.4 reveals 
the calculated solar radiation absorbed by the atmosphere and the 
earth's surface at the first of January. 

2.6.2 Longwave Radiation 

The calculation of the longwave radiative cooling makes use of a 

simple non-interactive scheme. The computations are performed by taking 

longwave flux differences across each layer. The upward and downward 

emissivities were estimated from values given by Stone and Manabe 

(1968). Those values were calculated using an 18 level radiation model. 

Average vertical profiles of temperature, water vapor, ozone and carbon 

dioxide were used for calculation of the longwave fluxes. 

The atmosphere in the present mooel is divided into two layers. 

The net longwave flux at the surface is given by 

L = £ B T (2.34) 

s s g 

The net longwave flux at 500 mb is 

L2 = £2 B V’ 

and the net longwave flux at the top of the model’s atmosphere is 

L 3 = B 12 '*. (2.36) 

B is the Stefen-Boltzman constant. £^, £2 and are the net emissivi- 
ties at the surface, 500 mb and the top of the atmosphere, respectively. 
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Table 2.1 

Parameters Used for Solar Radiation Calculations 

Parameter Clouds Atmospheric Surface 

Latitude c Albedo (r ) Albedo (r ) 

a s 


84.1 

0.35 

0.096 

0.8 

76.5 

0.41 

0.129 

0.8 

68.9 

0.48 

0.179 

0.8 

61.3 

0.54 

0.305 

0.4 

53.6 

0.56 

0.343 

0.3 

45.9 

0.54 

0.316 

0.2 

38.3 

0.45 

0.248 

0.2 

30.6 

0.37 

0.185 

0.18 

23. 

0.28 

0. 131 

0.15 

15.3 

0.29 

0.145 

0.14 

7.7 

0.32 

0.167 

0.14 

0. 

0.38 

0.207 

0.14 

-7.7 

0.36 

0.193 

0.12 

-15.3 

0.35 

0.183 

0.1 

-23.0 

0.34 

0.166 

0.1 

-30.6 

0.36 

0.179 

0.1 

-38.3 

0.42 

0.227 

0.1 

-45.9 

0.51 

0.293 

0.1 

-53.6 

0.60 

0.377 

0.5 

-61.3 

0.62 

0.369 

0.5 

-68.9 

0.55 

0.6 

0.6 

-76.5 

0.47 

0.6 

0.6 

-84.1 

0.40 

0.6 

0.6 
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0 

Latitude 


the earth tfuii iXw^ af»5phere (dotted U„o) and 
tne eartn Uuil line) for the first of January. 





39 


ORIGINAL PAGE 19 
OF POOR QUALITY 


Based on the Stone-Manabe model, the net emissivities were assigned 

the values 

£ = 0.26, 
s 

= 0.74, 
and = 0.11. 

Thus, the longwave heating rate of the lower layer is 

'■ 2 >' 

and the longwave heating rate of the upper layer is 

(QL )3 = Ip (L 2 - L 3 ). (2.38) 

2.7 Surface Temperature 

The bulk formulae (2.24) and (2.25) require the surface* temperature 
T^. Over land and ice surfaces the temperature is calculated from the 
thermal energy equilibrium on the surface, assuming negligible heat 
capacity of the earth (Holloway and Manabe, 1971). The net solar radia- 
tion absorbed by the earth’s surface (S^), the sensible heat exchange 
with the lower layer (Qg)> the longwave emission of the earth’s surface 
and the latent heat due to evaporation from the surface (Q^) ^te in 
balance; 

S = e BT ^ + Q + Q + I(T - 271.2°). (2.39) 

s s g g 

The last term on the right side represents the conduction of heat from 
unfrozen water below sea ice in the polar latitudes of the winter hemi- 
sphere. Assuming a thermal conductivity of ice, T^ =2.1 J m ^ ^ 

sec \ a temperature of the underlying water of 271. 2”K, and an ice 

2 -1 

depth, d = 2 m, then the constant I = ~ 1^05 W/m . This term 

is needed to prevent unrealistically cold temperatures in the Northern 
Hemisphere polar regions during winter. 
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Equation (2.39) is solved over land and sea ice areas using 
Newton’s iteration method (Appendix I). 

Over the oceans the sea surface temperatures are specified using, 
the mean January climatological values. Since the variation of sea 
surface temperatures between late winter and early spring is small, we 
keep these temperatures unchanged during the three month simulation 
(January, February and March). 

2.8 Large Scale Precipitation and Convective Adjustment 

The model has a moisture content in the lower layer (level 1) only. 

The procedure for large scale precipitation and the convective adjust^ 

ment starts after completing each time step of integration. The mixing 

ratio harmonics are transformed to physical space, and then the model's 

atmosphere is examined for super saturation at each grid point. 

If q(T-) < yq (T-), then no precipitation or convective adjustment 
X si 

takes place. The parameter y represents a specified critical relative 

humidity (y = 0.85 in this study), is the temperature at any grid 

point in level 1, and q and q are the mixing ratio and the saturation 

s 

mixing ratio, respectively. 

On the other hand, if q(T.) > yq (T,), condensation occurs with the 
associated latent heat release. The temperature Tj will be agumented by 
an increment AT, such that 

AT = ^ (q(Tj) - q^(T + AT)), (2.40) 

P 

where q' is the new saturation mixing ratio at the temperature T AT, 
s 


(2.41) 
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Using the Clausius-Clapeyron equation » (2.41) takes the form 




Lq 

Ft^ 

V 


AT 


(2.42) 


where is the water vapor gas constant" . The rate of condensation 
(precipitation) per unit mass, P^, is given by 


P = 

c 

where At 


P = 
c 


(q - q^Mt, 


is the time step of integration. 


q-Yq_ yL2 

At / c R Ti? ' 

p V 1 


Using (2.40) and (2,42) 
(2.43) 


It is clear that a relevant form of (2.40) is 


AT = — P At. (2.44) 

c c 
p 

After the release of latent heat in the lower layer as a result of 

the condensation of water vapor, the atmosphere is tested to see if 

convective adjustment is required. Convection is assumed to develop if 

the atmosphere is unstable relative to the moist adiabatic lapse rate 

r . If the temperature lapse rate is greater than P , then the tempera- 
s s 

ture of the two levels are adjusted to stabilize the model’s atmosphere 
by cooling the lower layer and warming the upper layer, with the verti- 
cally averaged temperature conserved. The new lapse rate is the same as 

r . 

s 

2.9 Model Equations 

The system of equations (2.1 - 2.6) together with the vertical 

boundary conditions discussed before (Section 2.3) are expanded at 
levels 1 and 3 of the model. Wc use the average values 4^ = (4^i + 4^s)/2 
and 0 = (0| + 03)/2 by adding the two systems of equations at the two 
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levels. The difference equations for l = (413 - 4 <i )/2 and o = (63 - 
6 i )/2 are calculated by subtracting the equations in both layers. For 


consistency of notation we introduce Xj such that X - Xj* 
system (2.1 - 2.6) will reduce to the following: 

Thus , the 

|^(V24() = -J(4,,v 2^+£) . J(x,v2i) - %V-fVXg+Fjj+F^, 

(2.45) 

|^(V2i) = -J(i|<,V2t) - J(l,V24)+f) + V.fVx -%V*fVXg+Fj^+F^, 

(2.46) 

36 

gr = -J(4<,6)-J(t,a)+V*(aVx )-^(Vx^ -VS+Vx ♦Va+3oV2x )+Q+W. , 
^ ^ s s s n 

(2.47) 

3o ^ ^ 

5T = -J(4<,a)-J(t,e)+Vx *V0-^(Vx ‘Ve+Vx *Va-oV2x )+Q+W. , 

(2.48) 

= -V*((KXV(||J-T)+Vx)q) + E - + (Sj^)j, 

(2.49) 

bCpV20 = V-(fVT), 

(2.50) 

and 


^^Xg = -2Vg ♦ Vn, 

(2.51) 

K K 

where b=l[[|] - [t] 1 = .124. 



This is a system of seven equations with seven unknowns 

<l», t, 0, a, Xi, Xg, q. 

All the budget values Fj^, F^^, F^, F^, W^, S^, E, P^, Q and Q are 

either specified or calculated, as discussed previously. 

^h ^ ^ = ’r ^ ( 2 . 52 ) 

Fjj + F^ = ^ - 2k^ V2x + ^ VH) (2.53) 

Wh = k^V2e 


( 2 . 54 ) 
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(2.55) 


= k,V^q 


(2.56) 


Q is the vertically averaged diabatic heating per unit mass, and Q 
is the difference in the diabatic heating between the two levels. 


‘liy vR 

’ll <■“> 

The surface stream function is extrapolated from the values at the 
two levels assuming that tjj is linear with height (Salmon and 
Hendershott, 1976). 

»|> = - 1.6t (2.59) 

s 

The above system represents the structure of the model. This system 
will be transformed to the spectral '"^ace using the spherical harmonics 
as basis functions. 


(2.57) 


2 . 10 Spectral and Finite Difference Methods 

The conventional spectral method is Galerkin's method based on ex- 
panding the different variables with a truncated series of surface 
spherical harmonics. The method is used for the numerical integration 
of the hydrodynamical equations. Two types of expansion are often used, 
the triangular and rhomboidal truncations. The advantages of the spec- 
tral method over the usual finite difference methods are summarized as 
follows (Machenhauer , 1974); 

1) The non-linear terms are alias free, which prohibits the exis- 
tence of the non-linear instability described by Phillips (1959). 
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2) Quadratic area integral invariants like the kinetic energy and 
enthalpy also are invariant for the truncated system, since the 
error fields are orthogonal to the variables. 

3) Linear terms are computed without any truncation error. 

4) No special treatment is required for dealing with the polar 
region when using the vorticity and divergence fields. By con* 
trast, in the finite difference method the horizontal wind com- 
ponents are discontinuous at the pole. 

5) The friction term of the finite difference methods is necessary 
to prevent aliasing instability. It also is necessary for the 
removal of energy from the shortwave end of the spectrum. When 
using the spectral method, it also is important to prevent blocking 
of energy at the highest wave numbers retained, but in this case 
the purpose is only a simulation of the effect of the small scales 
not retained in the representation. 

A study by Hoskins and Sim^cns (1974) compared finite difference 
and spectral models. The study showed that no one method has a super- 
iority in all respects. In comparison with the finite difference model, 
the spectral model gave much improved solutions for the amplitudes and 
phases of the predicted waves. On the other hand, the finite difference 
model gave a more accurate representation of the frontal systems. 

It is of interest to compare the two types of truncation mentioned 
before, namely the rhomboidal and triangular. For the same zonal wave 
number truncation, the triangular representation has fewer degrees of 
freedom than the rhomboidal and hence requires less computing time. If 
we retain the same degrees of freedom in both the triangular and rhom- 
boidal truncations, the former will be more appropriate for mean zonal 
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fields than the latter. At the same time the rhoaboidal truncation 
could introduce higher wave numbers, namely the eddies. The same study 
by Hoskins and Simmons (1974) did not give a definite conclusion con- 
cerning the comparison between rhomboidal and triangular truncation. In 
some experiments the rhomboidal truncation gave a more accurate approxi- 
mation to the solution than the triangular truncation. In other experi- 
ments the triangular truncation gave a more efficient description of 
Rossby wave instability. 

In this study we used the rhomboidal truncation since it gives a 
comparable resolution in both horizontal directions. 


2.11 Spectral Method 

The dependent variables X> X i <1 are expanded in trun- 

s 

cated series of the form 


M |m|+J 

X(M,X) =21 (|j,\) 

m=-M n=|ml 


n n 


(2,60) 


where X is any variable being studied, are harmonic coefficients, X 
is longitude, |j is the sine of latitude, m is the zonal wavenumber, n is 
the degree of a :>pherical harmonic component , | n- m| is a meridional 
wavenumber in the sense that there are|n- m| zero crossings of yJJ 
between equator and pole, M is the highest zonal wave number retained in 
the truncated series, and J is the highest value of(n- m| retained in the 
truncated series. are spherical harmonic functions defined by 


y“ = P™ (m) e 
n n 


imX 


(2.61) 


P™ are the Associated Legendre functions of the first kind 


."•r.n - / ill'll) (n~ n O ! 

^ 4/T ' (n+ ni)! 


,n+ m 

u 

, n^ m 
dM 


(M^-i) 


n 


(2.62) 
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A spherical harmonic coefficient is defined by 

- 2n +1 * 

X® = i- f f X y" 

*n . J J ^ dMd\ 

0 -1 

where Y^* is the complex conjugate of Y^. 


are orthogonal over 



the surface of the sphere, 

1 for (mj,nj) = (m,n) 

0 for (fflj.nj) / (m,n) ’ 


i.e. 


and are eigenfunctions of the Laplacian operator 


(2.63) 


(2.64) 



(2.65) 


where a is the radius of the sphere. The coefficients for negative and 
positive values of m are related in the following way; 

X"® = (-1)® X®*. 
n n 

Nonlinear terms are transformed from grid point space to spectral space 
using the full transform method (Machenhauer and Rasmussen, 1972; 
Orszag, 1970). The method is computationally highly efficient relative 
to the interaction coefficient method for J > 9. 

The procedure for calculating the spectral coefficients of the 
nonlinear terms using the full transform method is as follows: 

1) Calculate the nonlinear terms at each grid point in physical 
space. 

2) Transform to the Fourier space at each Gau^'8ian latitude, using 
fast Fourier transform routines. 

3) Transform to the spectral space using the Gaussian quadrature 


formula . 
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Highly nonlinear terms, like diabatic heating terms cause problems 
in finding their spectral transforms. This problem is resolved by using 
the full transform method. They are calculated in physical space, then 
added to the nonlinear dynamic terms, and the whole sum is transformed 
to spectral space. 

To guarantee an alias-free solution, there are two conditions that 
must be fulfilled (Machenhauer and Rasmussen, 1972). These conditions 
specify the minimum number of zonal grid points, N^, and the minimum 
number of Gaussian latitudes, I^, on the sphere: 

N > 3 M + 1 
g 

I > M + 3/2 J. 
s ' 

In case of the rhomboidal truncation (M = J) used here, the latter 

condition is 

I > 5/2 M. 
s 

For the simulation with wavenumber 9, N =32 and I =23. 

8 s 

To transform the system (2.45 • 2.51) to its spectral form, each 
variable is expanded using (2.60), Tae resulting equations are mul- 
tiplied by and integration of both sides is performed using 

equations (2.63 - 2.65). The nonlinear terms are calculated using the 
transform method mentioned before. 

The system of equations in its spectral form is given by 

^n n(n+l) n(n+l) 


^ .p+2 ,,m m ^ n-J ,jn in 

" <,ST “n.l u*l * “n „-l 


( 2 . 66 ) 


s , m ,, ( n+ 1 ) I in . h , m 

(j) - K, n 2 + -r- tj) 
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^ ,n>2 m ^ n-1 ^m m ^ 

" ‘s?f “ntl X, „*1 * -„- ”n X» n-l> 


u. s I n> ou ,™ u (n+1) ,m ^ n ,m 
+ ;p- * 2 k , T - k. n ' ■ V X + - 9 - X 

2 ^s n d n h a^ n a^ n 


(2.67) 


2fi (2l^ D™ y™ + D'” y™ ) 

^n+1 n+1 \+l n n \-r’ 


e“ = {-j(«i>,e) - J(x,o) - Jj(Vx -ve + vx -^o + sov^x )}“ 

n s s s n 


. „ 9* , (V(a7x));; * q;; . 


( 2 . 68 ) 


a™ = {- J(«|.,o) - J(x,6) - JiO'Xg-Ve + Vx^-Vo - oV^Xg)}" 


+ (Vx-ve)” - n k o” + Q® 

n h n n 


(2.69) 


q“ = - (V-((kxV(4,-x) + Vx)q))” - n q” + ';k-P^)“, (2.70) 


b c e" = 2Qi^ d” , xr, + ^ d” x" ,), 
p n n-^1 n^l n+1 n n n-1 


(2.71) 


4» “ = D»® - 1.6 x”, 

s n n n 


's n n(ni 


(v^-vn)^. 


where i = V"l- The spectral transform of terms of the form V‘(fVx) or 
^‘(f^X) is shown in Appendix (II). 


49 


It must be noted that by solving (2.67), (2.68) and (2.71) we can 
obtain an equation for x* The equations are simplified and solved as a 
system of tridiagonal matrices to find the spectral coefficients of x 
that satisfy the linear balance approximation. The simplification is 
needed to treat the term (V*oVx) in (2.68). To do this, we split a into 
its global average [a] , and the deviation from this average a* , 

a = [o] + a* . 

Then 

V-(oVx) = la)V2x + V-(a'Vx). 

The first term on the right side of the above equation is of a larger 
order of magnitude and is added to the other unknown terms, which in- 
clude X- The smaller, second term, is considered as a known parameter 
and calculated using the values of x at the previous time step. The 
method is found to be stable. It significantly reduces the number of 
calculations at this stage. 

2.12 Energetics of the Model 

The two layer model discussed here conserves the sum oi kinetic and 
available potential energy under reversil Ic adiabatic processes (Lorenz, 
196)). If one introduces the topographical forcing as a lover boundary 
ver:ical velocity, it is hard to verify the energy conservation (Burger 
and Riphagen, 1979). It is only the very simple lower boundary condi- 
tion = 0 (used by Lorenz) at p = 1000 mb that guarantees an energy- 
conserving system. 

The kinetic and available potential ener'^ies, K£ and AP, respec- 
tively, are expressed in the forms 
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and 


2b c ^ 
AP = ^ 


((e'llMoilil 


(2.75) 


^ la] + [a2 + ( 0')2 + (o’)2j^ 

The square brackets ( ] indicate the global area average and the 
indicate the deviation from that average. 

In spectral space the kinetic and available potential energy within 
a spherical harmonic mode are given by 

(KE)® = ^ ( |||“ 2 2) n(n+l) (2-6 ), for m > 0 (2.76) 

"^n ga^ ^ ^n n ^ om " - ' ' 

and 


(AP) 


m 

n 


2b c Ap ( 0® 2 + o® 2) (2-6 ) 

p ^ n n om 


S 


+ {II 2 + 



(2.77) 


rs 

for n ^ 0 m> 0, where 6 =1 and 6 = lor m > 0. 

— ^ oo om 

2.13 Initial Conditions and Time Integrations 

The model integration starts from a hypothetical, horizontally 

isothermal, atmosphere at rest with a moist adiabatic lapse rate. The 

model runs for sixty days assuming perpetual solar forcing (first of 

January) . This initialization procedure is used in order to reach a 

statistically steady state. After that the solar declination is changed 

daily to simulate the climates of January (days 61-90), February (days 

91-120), and March (days 121-150). These runs are considered as control 

runs for the comparable periods within the experiments. 

The time difference method used is the centered (leap-frog) scheme. 

To avoid the growth of unnecessary computational modes, a time smoother 

was used on the prognostic variables (Asselin, 1972) at every time step. 
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The diffusion terms were calculated using values at the previous time 
step to ensure computational stability. The time step used was 2 hours. 

2 , 14 The Quasi -Geost t ophi c Versus Primit ivt'_ 

The dominant energies of the ultra-longwaves (zonal wave numbers 
1-4) and longwaves (wave numbers 3-^) are quasi-geostrophic. The quasi- 
geostrophic approximation assumes near-bal.mce between the Coriolis 
force and the pressure gradient force. This balance between the two 
forces enables us to ignore the time variation ol the horizontal diver- 
gence among other terms in the divergence equation. The magnitude of 
the velocity potential (\) is about an order of magnitude less than that 
of the stream function (4’^- This alst> was notioed in results of the 
present model coiuenung the two variables. The relatively high fre- 
quency gravity waves art' iiltt'it'd. The filtering of the gravity waves 
gives the advantage of using a large time step when the prognostic 
equ.iiions are numerically uilegiated. Using the primitive equations 

requires a much smaller time stop to satisfy the stability conditions, 
unless the semi ~ imp I i c 1 1 scheme i * ust'd. Kveii with Itial proceduie, the 
time step is less than tin one used in liht'red models of the same 
resolution. The quas i ~gt'i>sl roph u modt'is have' bc't'u ust'd iov prediction 
of middle latitude c i rcu I ,it i ons wheit' llie u 1 1 i a- 1 ongwaves and longwaves 
are well defined by llu' balaiuo between the two forces mentioned above. 

2.14,1 Scale Analysis 

The conditions undt'r whicli the primitive' equations could be reduced 
to tlie quas i -ge'ost i opin c system can bt' summaii/eil as (t'.g., Phillips, 

1) H 1 

2 ) HU H-K V 1, 

i 


uib.n 
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where R and R. are the Rossby and Richardson numbers, respectively, and 
0 1 

BG is the Burger number. These nondimens ional numbers are defined as 
follows : 


R 


0 


_ U_ 
■ fL’ 



BG 




In these equations, U is the horizontal velocity; L is the horizontal 
scale; H is the vertical scale; Q is the angular velocity of the earth; 
f ;s the Coriolis parameter; g is the acceleration of gravity, and (^) 
refers to the vertical density stratification. The second condition 
impliss that R^ = -^r which cannot be fulfilled in the tropics since R^ 
ha: u typical value of nearly hO for atmospheric long and ultra long* 
wa' es . 

In this model the variability of f with latitude is retained. Due 
to the larger values of R^ in the equatorial regions, the quasi- 
geostrophic model cannot detect many features of the tropical circula- 
tion. 

2. A. 2 Effect of Filtering the Gravity Waves 

Besides the question of validity of the geostrophic approximation 
in tropical regions, it is also necessary to discuss the effect of 
filtering gravity waves on the tropical circulation. The mo't important 
type is the Kelvin wave. In the longitude-height plane a Kelvin wave is 
simply an ordinary internal gravity wave which propagates downward and 
eastward in the tropical stratosphere. Kelvin waves may provide a 
we:.terly momentum source for the semiannual oscillation. There is no 
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direct observational evidence linking the stratospheric Kelvin waves to 
any weather disturbances in the troposphere (Holton, 1976). There are 
some theoretical studies, using linear primitive equation models, chat 
tend to explain the east-west tropical Walker Circulation in terms of 
propagation of Kelvin waves (e,g., Gill, 1980). In all such studies, 
the Rayleigh friction in the troposphere was assumed to be very high as 
compared to the observed values. 

2.14.3 Validity of the Model for the Proposed Study 

Although a nonlinear quasi-geostrophic model will not represent the 
tropical circulation as accurately as a primitive-equation model, the 
former can be used efficiently as a general circulation model to study 
the response of the midlatitude atmosphere to lower boundary forcing in 
the tropics. A nonlinear quasi-geostrophic model will be more appli- 
cable than a linear primitive equation model for such studies. This 
statement is deemed true since the atmosphere itself is highly non- 
linear. Also, for a quasi-geostrophic model the fulfillment of the 
energy constraints and the large time step used are both desirable in 
long term integrations. Semi-implicit schemes give the advantage of 
using large time steps in the case of the primitive equation models. 
However, they still use smaller time steps than those used by a quasi- 
geostrophic model of the same resolution. The planetary Rossby waves, 
which are the major modes representing motions away from the equator, 
are quasi-geostrophic in character. At the same time, it is not yet 
clear how the Kelvin waves affect, or interact with, the lidlatitude 
planetary waves. In other words, the response of the midlatitude atmo- 
sphere to thermal heating in the tropics may not be affected by the 
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propagation of the Kelvin waves there. In the tropics the local re- 
sponse to the heating seems to be a direct type of circulation. In 
either case, many features of the tropical circulation could not be 
represented accurately with a two level model. 

Based on the previous discussion, the present model appears to be 
suitable and efficient for studying the response of the large scale 
midlatitude atmosphere to thermal anomalies in the tropical regions. 


3. COMPARISON OF THE MODEL *S WINTER CLIMATE WITH OBSERVATIONS 


In order to gain an understanding of the limitations of the models 
it is important to compare the model’s climate with the observed large 
scale circulation. The simulated January climate is the 30 day average 
of days 61 through 90. February and March simulations are the inte- 
grations from day 91 to day 120 and from day 121 to day 150, respec- 
tively. In this period of integration from day 61 to day 150, the solar 
declination is changed daily as a function of the mean zenith angle 
mentioned in Chapter 2. These simulated monthly averages are considered 
as control runs for the proposed anomaly experiments. 

In this chapter we will compare the simulated January circulation 
with the corresponding climatological fields for the same month. The 
Northern Hemisphere observed fields were taken from tables given by Oort 
and Rasmusson (1971). For some fields the values from 10°S to the south 
pole are taken from the observed December-February average values 
(Newell ^ , 1972). 

When comparing the simulated and observed fields we must keep in 
mind the typ^ of approximations and simplifications assumed in con- 
structing the model and also the aim and nature of this research. The 
dynamical simplifications include use of the linear balance model in 
order to filti^r out the gravity waves present in the system of primitive 
equations. Tie low order of truncation (truncated at zonal wave no. 9) 
excludes the higher wave numbers, which certainly affects the eddy sta- 
tistics of the model. 'Hie vertical resolution of two levels is the 
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minimum required to allow baroclinic development. The physical simpli- 
fications include using zonal average values for the radiative forcing. 
For example^ we assume a January climatological average for the surface 
albedo, atmospheric albedo, and cloud amount. The moisture budget was 
considered in the lower layer only. Such simplifications have their 
effects on the simulated climate of the model. Since we are mainly 
concerned with the response of the planetary and long waves in middle 
and high latitudes to different locations of the thermal anomalies, the 
model appears to be adequate The analysis of the experiments was done 
using the differences between the anomaly cases and the control run. 
Therefore the absolute values of the different fields are of secondary 
concern. 

3. 1 Average Zonal Values and Variances Resulting from Stationary and 
Transient Eddies 

In this section we will consider the 30 day averages, representing 
January, for different fields. The zonal averages and the variances 
resulting from stationary and transient eddies are compared with obser- 
vations. In a two level model, the tropospirric average for linear 
parameters is represented by the values on the intermediate level (500 
mb in the present model). Nonlinear parameters, like variances and 
transports, are calculated at each level; then are vertically averaged 
to represent the simulated tropospheric values. 

To clarify the analysis, we define the following statistical rela- 
tions for a certain variable X (Oort and Rasmusson, 1971): 

The time average 

- 1 ^2 

x = lj ;^xdt 


(3.1) 
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Departure from the time average 
X' = X - X 
Zonal Average 

1 ^ 

[X] = ^ / Xd\ 

0 

Departure from zonal average 
X* = X - [X] 

The variance of X resulting from transient eddies 


(3.2) 


(3.3) 


(3.4) 


TX = (3.5) 

The variance of X resulting from stationary eddies 

SX = [X^'2] (3,6) 

The calculation of the transient variances gives an indication of 
the development of disturbances in the model (i.e. the weather associ- 
ated with the month of simulation). Naturally this variance will be 
more dominant in the middle and eytratropical latitudes of the winter 
hemisphere. The stationary eddy component reflects the effects of the 
permanent forcing. This results mainly from the diabatic heating dif- 
ferences due to land and sea distributions and also the topographical 
gradients. 

3.1.1 Zonal Wind Component 

In the Northern Hemisphere the simulated 500 mb mean zonal wind is 
stronger than the observed values (Fig. ;K1). The maximum value appears 
at 30®N, in agreement with observations, but it is 6 m sec ^ more than 
the observed maximum. In tropical latitudes the simulated belt of east- 
erlies is shifted to the Southern Hemisphere, with an amplitude about 4 
m sec ^ stronger than the observed. The simulation of the easterlies in 
the winter hemisphere seems to be a problem when using balance or quasi- 
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geostrophic dynamics (e.g. Wielicki and Hendershott, 1979; Kikuchi, 
1969, 1979). In the Southern Hemisphere the simulated maximum wester- 
lies are less than those revealed by observations. The characteristics 
of the simulated pattern in both hemispheres may be a result of the 
simulated temperature fields and the associated thermal wind components. 

In general, the variance in the horizontal wind field is a measure 
of the eddy kinetic energy in the atrrosph(?re , and, therefore, it is one 
of the basic parameters in the description of the climate. The sim- 
ulated, vertically averaged variance of zonal wind resulting from the 
transient eddies has less amplitude than the observed pattern (Fig. 
3.2). In the Northern Hemisphere, both the simulated and observed 
profiles have the same latitude of maximum variability (about 50°N). In 
the Southern Hemisphere (summer hemisphere) the simulated variance is 
generally small. The variance resulting from the stationary eddies 
(Fig. 3.3), has an amplitude which is comparable with observations, 
although its maximum is shifted about 10® equatorward within the 
Northern Hemisphere. Although the observational data are not available, 
the simulated stationary variances in the southern hemisphere seem to be 
in the correct position. Also they have less amplitude than the simu- 
lated Northern Hemisphere values. 

3.1.2 Meridional Wind Component 

The mean meridional circulation plays a central role in the opera- 
tion of the general circulation of the atmosphere. Unfortunately, as 
stated by Oort and Rasmussen (1971), it has been a very difficult task 
to obtain a reliable estimate of the intensity of the mean meridional 
circulation. This results from the fact that the zonally averaged mer- 
idional wind component must be estimated from the divergent part of the 
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wind field, which is generally of a lower order cf magnitude than the 
rotational part in the middle and high latitudes. At the same time, in 
the tropical region, the zonally averaged velocity potential is hard to 
detect in the observed wind field?. The ^bseived maximum 500 mb mean 
meridional wind component is about four tii*^es larger than the simulated 
value (Fig. 3.4), but still less than 1 m sec In the Northern 
Hemisphere, the simulated maximum southward mass transport is at 38.5®N, 
whi*'h agrees with observations. In the Southern Hemisphere the simu- 
lated pattern changes from southerlies to northerlies at 20®S, whereas 
observations show the change at 30®S. The small values of the simulated 
meridional wind component are partially a result of zonal average 
forcing used in the model, which underestimates the actual east-west 
differences . 

The calculated variance due to transient eddies (Fig. 3.5) has the 
same latitudinal position of the maximum as the observed one, but with 
less amplitude. Both, the position and amplitude of the variance re- 
sulting from stationary eddies (Fig. 3.6) show a better agreement with 
observation. The positions of the simulated values in the Southern 
hemisphere also seems to be in reasonable agreement. 

It is of interest to note, the close similarity of the patterns of 

(u’^1 and (v*^) for both the simulated and observed transient eddies. 
This indicates that the kinetic energy of the developing disturbances 
seems to be almost equally distributed among the zonal and meridional 
wind components At the same time, both the simulated and observed 
fields do not show such similarity in the case of variances resulting 
from stationary eddies. In the Northern Hemisphere middle and extra- 
tropical latitudes, both the model and observations reveal that the 
transient components are larger than the standing components. 
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and simulated (dotted line) zonally averaged variance 

2 “ ^ 

of zonal wind (m sec resulting from stationary 
eddies for January. 



The zonal average of the observed (full line) and 

simulated (dotted line) meridional wind (m sec ) at 
500 mb for January. 
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The tropospheric average of thf* observed (full lire) 
and simulated (dotted line) zonally averaged variance 

2 •‘2 . 

of meridional wind (m sec ) resulting from transient 
eddies for January. 


Stmiiiatad 


Obsarvad 


The tropospheric average of the observed (full line) 
and simulated (dotted line) 2:onally averaged variance 

2 -2 

of meridional wind (m sec ) resulting from stationary 
eddies for January. 
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The simulated variance of horizontal wind resulting from transient 
eddies is smaller than that revealed by observations (nearly about 1 to 
6). A part of this deficiency is due to the north-south horizontal 
resolution used in the model (Manabe et al. , 1965). 

3.1,3 Vertical Pressv.re Velocity 

The observed values were calculated assuming a zero lower boundary 
pressure velocity. In other words, the effect of topography was not 
taken into account for the observed fields. In the Northern Hemisphere 
the major areas of subsidence appear to be in the subtropics for both 
the simulated and observed fields (Fig. 3.7). The simulated rising 
motion in the tropical Southern Hemisphere associated with the active 
convection there is broadened, but with lower magnitude, than shown by 
observations. However, the main differences between the observed and 
simulated zonal average vertical velocity are apparent in the extra- 
tropical latitudes of both hemispheres. In the winter hemisphere the 
latitude of maximum vertical motion seems to be shifted poleward about 
12® from the observed maximum. 

The geographic map for the simulated vertical motion (Fig. 3.8) 
shows that the model correctly simulates the position of the maximum 
rising motion over the eastern Pacific and Atlantic Oceans due to the 
effects of the Rocky Mountains and Alps, respectively. The Himalayas 
cause a rather noisy pattern over eastern Eurasia and the western 
Pacific. In the Southern Hemisphere the effect of the Andes also is 
clear over the South American continent. One should also note, the 
upward motion over the tropical western Pacific and the sinking motions 
over the eastern Pacific adjacent "o the South American continent. This 
mi.cion pattern may he interpreted as the vertical branch of a Walker 
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The zonal average of the observed (full line) and 
simulated (dotted line) vertical pressure velocity 
-4 -1 

(10 mb sec ) at 500 mb for January (negative values 
upward, positive values downward). 
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3.8 The simulated vertical pressure velocity at 500 mb 

(10~^ mb sec”S for January (negative values upward, 
positive values downward). 
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type circulation, although the Kelvin waves were filtered from the 
model’s equations. 

The zonal average variances resulting from stationary eddies (Fig. 
3.9) show that the maximum variability of both the observed and sim- 
ulated patterns is between 20®N and 40®N. The simulated pattern does 
not show the secondary maxima near the equator and 70®N. In the 
Southern Hemisphere the simulated values give a small maximiun between 
20^3 and S0®S. 

3.1.4 Fiee Air Temperoture 

Tvhe zonal average of temperatures at 500 mb is shown in Fig. 3.10. 
The January simulated field is colder than the observed. In addition 
the temperature difference between the equator and 75®N exceeds the 
observed difference by about 10®C. Both patterns are relatively flat 
between lO^S and 40®N. However, north of 40°N the Simula teu temperature 
gradient is steeper than revealed by the observations. In the Southern 
Hemisphere we cannot rely on the observed (December - February) averages 
since January may differ from the three-month average. In general, the 
simulated temperature is warmer than the observed between 35®S and 7S®S. 

The geographic' maps of the simulated January temperatures for the 
750 mb and 500 mb levels (Fig. 3.11a,b) reveal a reasonable simulation 
of th^se fields. In particular, note the configuration of the thermal 
ridge to the west of the North America Continent and Eurasia, and also 
the thermal troughs along the east coasts of these continents. In the 
Southern Hemisphere the warm centers over South America and Australia 
are correctly simulated. 

Statistically speaking, the variance of temperature is a measure of 
the deviation of the individual temperature distributions from the zonal 
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3.10 The zonal average of the observed (full line) and 
simulated (dotted line) temperature (®C) at 500 mb for 
January. 
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mean value. From a physical point of view, the variance of temperature 
is a measure of the local contribution to the eddy available potential 
energy (Lore z, 1955). The simulated tropospheric average of tempera- 
ture variance resulting from transient eddies (Fig. 3.12) shows a maxi- 
mum at 65®N, which agrees with the corresponding observed pattern. The 
amplitude of the observed maximum is about two times larger than the 
simulated one The maximum variance occurs at this latitude as a result 
of the high variability of temperature over the Eurasian and North 
American continents. Over these continents large variation in tempera- 
ture can occur due to alternative advection of polar and tropical air 
masses. In the Southern Hemisphere the simulated variances are smaller 
than those over the Northern Hemisphere, with a maximum value near 50®S. 
The variance of the simulated temperatures resulting from stationary 
eddies (Fig. 3.13) shows a good agreement with observations for both the 
location and amplitude. Both reveal a winter maximum between 50^N and 
60®N. This maximum occurs as a result of the contrast between warm 
oceans and cold continents. 

3.1.5 Surface Air Temperature 

The simulated (1000 mb) and observed January surface air tempera- 
tures are shown in Fig. 14a, b. The large scale temperature distribution 
is well simulated by the model. The reservoirs of relatively cold air 
over the continents in the Northern (winter) Hemisphere and the rela- 
tively varmer air over the Southern (summer) Hemisphere continents are 
reasonably well projected. The 0®C isotherm is nearly in its observed 
mean position over the continents ir the simulated field, whereas over 
the oceans it is shifted southward from its observed position. This 
indicates that the simulated amplitude is less than the observed, as a 
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The tropospheric average of the observed (full line) 
and simulated (dotted line) zonally averaged variance 
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The tropospheric average of the observed (full line) 
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of temperature (^C^) resulting from stationary eddies 
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result of the zonal average radiative forcing used n the model. The 
cold pools over the eastern United States and Eurasia and the cold air 
mass associated with the Siberian high also are reasonably well simu- 
lated. In the Southern Hemisphere, the temperatures are somewhat colder 
than those revealed by observations. 

3.1.6 Geopotential Height 

The zonal average 500 mb geopotential height departures from the 
global average (Fig. 3.15) show the model’s tendency to predict lower 
pressure in the Northern Hemisphere and higher pressure in the Southern 
Hemisphere. The larger simulated pressure gradient in the Northern 
Hemisphere also is shown. The simulated January 750 mb and 500 mb 
geopotentials are compared with the observed 800 mb and 400 mb geopo- 
tential (Fig. 3.l6a,b; Fig. 3.17a,b). We can summarize the comparison 
as follows: 

1) The ridges on the western sides of North America and Eurasia 
are correctly simulated for both the "'SO mb and 500 mb fields. 

2) The ridge north of the Himalayas is simulated in the 750 mb 
map, as is the trough to the south, but the model calculates a 
higher geopotential gradient in the southern part, which leads to 
an over-estimate of the strength of the southwesterly flow there. 

3) The Aleutian low covers the central and eastern Pacific in the 
750 mb simulation. On the 500 mb map, the trough over the Pacific 
has a reasonable phase, but a smiller amplitude than that observed 
at 400 mb. 

4) In the Southern Hjmisphere the simulated 750 mb geopotential 
height shows the subtropical high pressure region over the eastern 
Pacific and western Atlantic. The observed low pressure area 
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The geopotential height for January. a) observed 
(lO^m) at 800 mb; b) simulated (ra) at 750 mb. 
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between 10**N and 10*’S in the Pacific is reasonably siaulated. The 
sinulated 500 nb geopotential shows a reasonable agreeaent with the 
observed 400 ab field. The latitude of the siaulated aaxiaua 
geopotential gradient is at about 50”S, in agreeaent with the 
observations . 

In the Northern Heaisphere the location of the variance of the 
siaulated geopotential heights resulting froa transient eddies agrees 
with the observations, but with a saaller aagnitude (Fig. 3.18). The 
observed aaxiaua, between latitudes SO^N and 75*’N, is about three tiaes 
larger than the simulated maximua. The same conclusion is true for the 
variance resulting from stationary eddies (Fig. ^.19), but the simulated 
amplitudes are about half the observed values. 

3.1.7 Sea level pressure 

Although the sea level pressure may be not the most important 
factor from the point of view of climate studies, its representation 
completes a summary of the general circulation. The model does not 
predict the surface pressure. A constant lower boundary of 1000 mb is 
assumed. However, we can still have an estimation of the sea surface 
pressure of the model using a diagnostic relation that takes into con- 
sideration the average geopotential, temperature and the static stabil- 
ity of the model's atmosphere. This relation is derived, using the 
hydrostatic relation, and assuming an atmosphere with constant lapse 
rate, 

Z, g/RF 

p = Ap (1 + r =^) , (3.7) 

s 

where P is tie temperature lapse rate, p is an analogue to the surface 

s 

pressure (mb), Ap = 500 mb, Z 2 is the 500 mb geopotential height, and T 2 
the 500 mb temperature (°K). 

c 
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3.18 The tropospheric average of the observed (full line) 
and simulated (dotted line) zonally averaged variance 
of geopotential height (m^) resulting fron transient 
eddies for January. 



3.19 The tropospheric average of the observed (full line) 
and simulated (dotted line) zonally averaged variance 
of geopotential height (m^) resulting from stationary 
eddies for January. 
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Fig. 3.20 represents the January sinulated sea level pressure 
deviations from a global avenge. ' In ’general we can see the tendency 
for the model to simulate a mean sea level pressure which is low in the 
winter and high in suoMr. This tendency was mentione^^ ' y Gates and 
S^hlesinger (1977), who used a two-level primitive equal Ion nodel with a 
sigma coordinate system. In the Northern Hemisphere the model correctly 
locates the Siberian high, the ridge over the western United States, and 
the Aleutian and Atlantic lows. Neither the Aleutian nor Atlantic lows 
spread over the entire ocean areas, but they do cover the central and 
eastern parts. In the Southern Hemisphere high pressure ridges appear 
over the subtropical South Pacific, the South Atlantic and che Indian 
Oceans. The relatively low pressure over the continents of South 
America and Australia also are simulated, while over South Africa the 
observed trough is shifted to the west. 

3.1.8 Mixing Ratio 

The moisture budget is calculated in the lower layer of the model 
only (750 mb level). Since the mixing ratio and the specific humidity 
are nearly the same and have the same properties, the simulated 750 mb 
mixing ratio for January is compared with observed specific humidities 
at 700 mb and 850 rot The zonally averaged simulated mixing ratio (Fig, 
3.21) shows the same tendency as the observed field at all latitudes, 
but is underestimated in the extratropical latitudes, where colder 
temperatures are simulated. 

The average tropospheric variance of the simulated and observed 
mixing ratio resulting from both the transient and stationary eddies are 
shown in Figs, 3,22 and 3,23, respectively. The simulated transient 
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3.21 The zonal average of the observed (700 mb, full line; 

850 mb, dashed line) and simulated (dotted line)roixing 
ratio (gm/kg) for January. 
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The tropospheric average of the observed (full line) 
and simulated (dotted line) zonally averaged variance 
of mixing ratio (gm/kg) ^ resulting from transient 
eddies for January. 
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The tropospheric average of the observed (full line) 
and simulated (dotted line) zonally averaged variance 
of mixing ratio (gm/kg)^ resulting from stationary 
eddies for January. 
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component has its maximum at IS^N^ abouc 5^ north of the observed maxi-- 
mum. Also, in the Southern Hemisphere, a secondary maximum is simulated 
at 40^S. The simulated stationary component has its maximum at 30^N, 
whereas the observed maximum is at 20®N. Both the simulated and ob- 
served fields reveal •'n area of maximum variability beU/een 10®N and 
40®N. There is a secondary maximum in the Southern Hemisphere near 
40®S. 

3.1.9 Precipitation and Evaporation 

Precipitation occurs in the model as a result of condensation due 
to large scale supersaturation. All condensed water is assumed to fall 
instantaneously to the surface as rain, without being re-evaporated or 
stored in the lower layer. Convection develops if the resulting lapse 
rate is larger than the moist adiabatic lapse rate, after the release of 
latent heat in the lower layer. 

The model’s simulation of precipitation seems to be more successful 
in the Northern Hemisphere than in the Southern Hemisphere (Fig. 3.24). 
In both hemispheres, more precipitation is simulated over the oceans 
than over the continents. According to observations, this is only true 
in the winter hemisphere (see for example Gates and Schlesinger, 1977). 
The zonal average of the simulated precipitation (Fig. 3.26) shows 
maxima on the equator, and at 30®N, 30®S, 70°N c'nd 80®S. The excessive 
^-ain in the Southern Hemispheric polar latitudes is considered to be in 
error. 

The evaporation from the surface to the model’s lower layer is 
better simulated than the precipitation field (Fig. 3.25). In both 
hemispheres the maximum evaporation is over the oceans, in agreement 
with observations (Gates and Schlesinger, 1977). The areas of maximum 


84 



ORIGINAL page 13 

nil^LlTY 



3.24 The simulated precipitation rate (cm/day) fcr January. 
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evaporation over the western Pacific, western Atlantic, Indian Ocean and 
South Pacific are reasonably well sinulated. One of the reasons that 
the simulated evaporation has a better match with o’ lervations than 
precipitation is that the model uses a similar bulk tormula as in cal- 
culating the evaporation from observed data. In parameterizing pre- 
cipitation, many important mechanisms were excluded, in addition to 
assuming that the moisture w^.s in one layer only. 

The difference between evaporation and precipitation represents the 
net moi..tur« source or sink to the model's atmosphere. The zonally 
averaged simulated sources (evaporation larger than precipitation) of 
moisture are in the subtropical latitudes of buth hemispheres between 
5°N and 28°N, and between lO^S and 25°S. Another source appears in the 
Southern Hemisphere between 45**S and 60°S. This source represents the 
effect of large amounts of evaporation from the ocean during summer. 
The sinks appear mainly in the equatorial belt between 10^’S and S*’N, 
55^N and the North Pole, and 25°S and 45**S. The large negative values 
near the South Pole, which result from the spurious amount of rain 
simulated by the model at this latitude, do not affect the results 
significantly due the small area of the polar cap as compared to other 
latitude belts 

3.2 Meridional Transports 

For the atmosphere-oce)>:i system, the zonally averaged radiation 
pattern shows an excess of incoming over outgoing radiation in low lati- 
tudes and a deficit in high latitudes. The radiational imbalance is 
maintained by strong poleward transports of energy in the atmosphere and 
oceans. The major components of energy transported are sensible heat, 


88 


ORIGINAL PAGE 18 


potential energy, latent heat and kinetic energy. In a certain lati- 
tudinal band, the time averages of these energy components usually 
balance. The convergence and divergence of the north-south flow com- 
ponent play an imporcant role in creating such a balance. We can define 
two mechanisms of transfer; namely, the transfer by the eddies (tran- 
sient or stationary) and the transfer by the mean meridional circula- 
tion. 

Using the same definitions as in section 3.1, we define the fol- 
lowing transport processes for a certain parameter X: 

northward transport by transient eddies = [v'X'] 

northward transport by stationary eddies = [^X'^] 

northward transport by mean meridional circulation = [v]''(X], 


where v is the meridional wind component and the superscript " denotes 
the deviaticn from the vertical average. 

The simulated transports for January of westerly momentum, sensible 
heat, potential energy and latent heat are discussed and coispared with 
the observed values. Kinetic energy transports will not be presented 
here. In general their influence is much le'js than that due to sensible 
heat, latent heat and potential energy (Starr, 1951). 


3.2.1 Transport of Westerly Momentum 

Both the observed and simulated meridional transports of westerly 
momentum by the transient eddies (Fig. 3.27) show a maximum near 35^N to 
45”N. This maximum is a result of the strong westerly subtropical jet 
stream. As in the case of the variances, the simulated amplitude is 
consid ‘rably smaller than the observed one. The transport of westerly 
mome..«.iiffl by stationary eddies (Fig. 3.28) shows a better agreement with 
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the observed pattern than the transient case. The maximum transport 
appears to be between 30®N and 40®N. The observed strong southward 
transport in the extratropical Northern Hemisphere latitudes is simu- 
lated as a small southward transport near the polar region. The model 
fails to Simulate the small southward momentum transport in the Northern 
Hemisphere tropical region. Jn the Southern Hemisphere there is a 
simulated maximum near 50®S. The observed transport of westerly momen- 
tum by the mean meridional circulation (Fig. 3.29) shows four main 
latitudinal zones with alternating signs. These are the southward 
transports in the Southern Hemisphere tropical region and iu the sub- 
tropical to midlatitude region, and the northward transport in low 
latitudes and the polar region. These four cells of transport are sim- 
ulated by the model, but in the area north of 30°N the model did not 
correctly simulate the location of the observed pattern. In the 
Southern Hemisphere the small observed southward transport in the tro- 
pical region is extended to reach the Southern Hemisphere subtropics in 
the simulated case. 


3.2.2 Transport of Sensible Heat 

In discussing this type of transport we must take into considera- 
tion that the sensible heat by itself is not a conservative quantity. 
Lament heat, potential energy and kinetic energy are needed to attain a 
total energy conservation. The latitudinal cross section of the simu- 
lated transports at 750 mb and 250 mb are compared with the observed 
cross sections at 700 mb, JOO mb and 300 mb. Observations show a maxi- 
m - jnsport by transient eddies (Fig. 3.30a) between 40°N to 60^N at 
ail lev s. The simulated pattern shows the maximum in the same lati- 
tudinal belt Fig. (3.31a), but with a slightly smaller magnitude. At 
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250 mb the simulated values show a small southward transport at most of 
the subtropical latitudes of the Northern Hemisphere* The simulated 
transport of sensible heat by standing eddies (Fig* 3.31b) has a reason- 
able agreement with the observed cross section (Fig. 3.30b), especially 
in the middle and extratropical latitudes of the Northern Hemisphere. 
In the Southern Hemispheric midlatitudes a simulated northward sensible 
heat transport appears at all levels. The transport by the mean meridi- 
onal circulation (Fig. 3.32) shows nearly the same differences as the 
transport of the westerly momentum. The simulated transport in the 
tropical region is extended to the subtropics in the Southern Hemisphere 
and to the midlatitudes in the Northern Hemisphere. In the extra- 
tropical latitudes of the Northern Hemisphere, the two patterns are not 
similar. One should emphasize the small values of this type of trans- 
port as compared to the correctly simulated transport by the eddies. 

3.2.3 Transport of Potential Energy 

The observed transport by transient eddies (Fig, 3.33a) shows a 
maximum southward transport in the middle and upper troposphere of the 
Northern Hemisphere between 30®N anr" 50°N. Also there is a secondary 
maximum of northward transport between 60®N and 80®N. These features 
are well duplicated by the model (Fig. 3.34a). At the lower level (750 
mb), the observed northward transport is stronger than the simulated 
transport. In the Southern Hemisphere the model simulates a southward 
transport in the lower layer and northward transport in the upper layer. 
The transports by stationary eddies show a similar agreement between 
observed i id simulated cross sections (Figs, 3.33b, 3.34b), particularly 
in the ^rthern Hemisphere subtropical and middle latitudes. In high 
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3.32 The tropospheric average of the observed (full line) 
and simulated (dotted line) zonally averaged northward 

transport of temperature (m sec"^ °C) by the mean 
meridional circulation for January* 
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Vertical cross section of the observed 2onally averaged 

2 

northward transport of geopotential height (m sec ) 
for January* a) by the transient eddies; b) by the 
stationary eddies. 
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The tropospheric average of the observed (full line) 
and sinmlted (dotted line) zonally averaged northward 

2 •! 

transport of geopotential height (m sec ) by the mean 
meridional circulation for January. 
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The tropospheric average of the observed (full line) 
and simulated (dotted line) zonally averaged northward 

transport of water vapor (m sec ^ gm/kg) for January, 
a) by the transient eddies; b) by the stationary 
eddies; c) by the mean meridional circulation 
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latitudes the observed upper-tropospheric northward transport is simu- 
lated, but with lesser magnitude, Th^ transport by the mean meridional 
circulation (Fig. 3.35) is similar to the transport of the sensible 
heat. The northward transport of potent* al energy in the tropical 
region is extended to the middle latitudes. A small southward transport 
in the extratropical latitudes and northward transport in polar lati- 
tudes are simulated. In the Southern Hemisphere, the model simulates a 
continuous southward transport by the meridional circulation south of 
35 »S. 

3. 2. A Transport of Latent Heat 

The transport of water vapor resembles that of latent heat. Al- 
though the hydrologic cycle used in the model is extremely simple, it is 
able to simulate the major observed features of latent heat transport. 
The simulated transport by transient eddies has a maximum between 20**N 
and 35"N (Fig. 3.36a), but with less amplitude than the observed pat- 
tern. The observed maximum transport by stationary eddies at 50°N (Fig. 
3.36b) is simulated at 30”. For both types of eddies, the model simu- 
lates a southward transport in the Northern Hemispheric tropical region, 
whereas the observations show this transport to be farther to the south. 
The simulated transport by the mean meridional circulation (Fig. 3.36c) 
also shows a reasonable agreement with observations, especially in the 
tropical and subtropical latitudes. 

3.3 S’ unmary 

. spite of its simple structure, the model could simulate many of 
the .aj..'r large scale features of the observed January circulation. We 
must tak-' into consideration that the observed statistics used for 
veiificat. iclude all the observed wave lengths. A fairer comparison 
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between the statistics of both sets would be made if ve were to truncate 
the observations at the same wave number used in tt iodel ( ive no. 9). 

From the above discussion it can be seen that the m< lei has two 
main deficiencies - the first in the precipitation field (near the poles 
and in the Southern Hemisphere) , and the second in the transient eddy 
motions. The simple hydrological cycle, with a crude precipitation 
scheme, is one factor in causing such deficiencies. The spectral trans- 
form method has its difficulty in simulating small quantities like the 
mixing ratio, particularly in the case of low ord'r truncation. Har- 
monic coefficients may yield negative values in the mixing ratio near 
the poles. This difficulty is treated by assuming a spurious source of 
evaporation to increase the mixing ratio to a i-±ro value. This false 
source of evaporation is balanced by an equal amount of precipitation. 
The small negative mixing ratios are set to zero. In general the effect 
of moisture on th(' general circulation is very small near the poles 
where such a procedure is necessary. 

The simulated values of the variability and transport caused by the 
transient eddies generally were underestimated by the model. This is 
expected from this type of simulation because of the limited number of 
layers representing the atmosphere, and, also, because of the low order 
of truncation assumed. Other points to be considered are the effect of 
diffusion of various scalar parameters and the time filtering scheme 
used, which were necessary for computational stability. At the same 
time the locations of such variances or meridional transports by the 
transient eddies usually were in reasonable agreement with observations, 
especially in extratropical latitudes, where such mechanisms are impor- 
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Similar problems were encountered using a more complicated model 
(Manabe ^ al. , 1965). The poleward transport of angular momentum by 
the large scale eddies in Manabe 's model was nearly half the amount 
revealed by observations. The horizontal resolution was found to be a 
major factor in improving these two deficiences. When the pole to 
equator grid intervals were increased from 20 to 40, the results were 
greatly improved, especially with regard to the vertical motion and 
precipitation. It is of interest that the dry version of Manabe* s model 
revealed greater transient transports than the moist version of the same 
model. The release of latent heat of condensation decreases the tem- 
perature gradient and the vertical wind shear in middle latitudes, 
decreasing the eddy kinetic energy in middle latitudes, and accordingly 
decreasing the northward fluxes of heat at.d momentum due to large scale 
eddies. 

The simulated variances especially in the temperature and geopoten- 
tial fields resulting from the stationary eddies show a reasonable 
agreement with observations in both, phase and amplitude. This is to be 
expected, since the ocean-land differences and the topographic forcing 
were included in the model. 

The transports by the mean meridional circulation were correctly 
simulated in the tropical and subtropical regions, whereas in the middle 
and extratropical latitudes, where the influence of such transports is 
generally small relative to that by the eddies, the simulated patterns 
do not reveal the same phase as the observ.^tions. 

Tn all cases the difference between the simulated amplitudes of 
di cerent fields between the Northern and Southern Hemispheres were 
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The previous discussion shows that the aiodel could te used for 
climate sensitivity <itudies. The main issue of the present study is to 
find the varial«ility 'n the quasi-stationary ultra-long waves and long 
waves as a result of thermal anomalies. These types of waves are well 
represented by both the dynamic structure and the horizontal resolution 
of the model. However, the model is not appropriate for the study of 
smaller-scale phenomena. 


ORIWNAL pa« « 


4. EX»’l':RIMi'.NT DESIGN AND STATISTICAL METHODS 
4 . 1 Inl roduct ion 

The experiments were designed to investigate the possible telecon- 
nections between local anomalies and the atmospheric circulation in 
subsequent months or seasons. As was discussed in Chapter 1, observa- 
tional studies, supported by many niunerical modelling results, suggested 
that tropical SSTAs in the Pacific ocean have significant effects on the 
atmospheric circulation in middle end extratropical latitudes, especial- 
ly in the Northern Hemisphere. Although many observational studies also 
suggested that middle latiluile SSTAs could have a significant effect on 
the atmosphere there, most of the numerical experiments failed to sup- 
port such a hypothesis. The negative phase of the Southern Osc'llation 
is defined as low surface pressure anomalies in the tropical east 
Pacific and high surface pressure anomalies in the Indonesian region. A 
certain SSTA pattern was found to be .associated with the negative phase 
of the Southern Oscillation. This pattern may be identified by the 
presence of warm water in the equatorial east Pacific (El Nino) and cold 
water in the central and western parts of the North Pacific. Horel and 
Walla..' (1981) suggested that this phase of the Southern Oscillation 
is a;;ociated with an anomalous mid-tropospheric pressure pattern in 
t,ne iddle and extratropical latitudes. The major features of this 
mid-tropuspheric pressure anomaly are high pressure in western Canada 
and low pressure over the southeastern United States and in the North 
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Pacific. This distribution was supported by the results of experioents 
using a linear model (Hoskins and Karoly, 1981). 

One phase of the experiments performed in this study is the inves- 
tigation of the effect of the SSTA in the Pacific Ocean using a more 
realistic nonlinear model with topographic forcing. The time changes of 
nonlinear interactions are important in middle and extratropical lati- 
tudes, where the advection of atmospheric variables is of a significant 
magnitude . 

Another phase of the experiments is to examine the response of the 
atmosphere to mid-tropospheric thermal anomalies . It is known that 
surface temperature has an effect on the atmosphere through the flux of 
moisture and sensible heat into the boundary layer. Host numerical 
models do not have a detailed structure of the boundary layer, and even 
those which have a boundary layer fail to represent many processes that 
take place there. The details of the complex mechanisms by which such 
fluxes affect the static stability of the lower troposphere are not 
fully understood. Because of crude parameterizations, the thermal 
effects caused by lower boundary anomalies may not resemble the actual 
tropospheric response. El Nino events, with warm water in the tropical 
eastern Pacific, are associated with cloud clusters across the width of 
the equatorial troposphere in the central Pacific region. If the SSTAs 
were inserted in a complex atmospheric model, they might not actually 
reproduce the same cloudiness structure, and, consequently, the same 
tropospheric thermal influence. Here, we assume that a middle tropo- 
spheric thermal forcing exists at certain locations. The anomalies are 
represented by propagating waves, which simulate the easterly waves in 
the tropical troposphere. 
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Exp e rimental D^si^n 

Starting from horizontally isothermal atmosphere initially at rest 
and with a constant lapse rate, the model was run for 60 days. During 
this period the solar declination of the sun was held constant in order 
to represent the first day of January. A statistically steady state is 
reached during these two months of integration. The simulated month of 
January is obtained by running the model from day 61 to day 90, the 
February simulation is the 30 day average from day 91 to day 120, and 
finally, the March simulation is the average of days 121 to 150. During 
the three month simulation, solar radiation is changed daily using a 
mean zenith angle. The simulation during these months has been done 
without including any anomalies. The sea surface temperatures are kept 
unchanged, representing the average observed values for January. The 
monthly simulations described above represent the control run of the 
model . 

The results of the anomaly experiments are analyzed and compared 
with the appropriate month of the control run. Each experiment is 
started from the initial state on the first of January (day 1). The 
anomalies are introduced gradually with an e-folding time of 5 days in 
January until they reach their full amplitude through the remaining days 
of January, February and March. An additional run is made with the SSTA 
turned off at the beginning of March. In this run the anomalies are 
gradually removed with the same enfolding time (5 days) with which they 
were initiated at the beginning of January. The differences between the 
March run with and without anomalies will give an estimate of the 
’’memory** of th atmosphere to different anomalies. 
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4.3 Statistical Significance of the Signal 

The purpose of any statistical significance testing is to provide 
an estimate of the probability that a given experimental resi It can be 
attributed to a specific forcing or perturbation applied ^ * the experi- 
ment. In general, any changes in the initial cond^'^ons or in the 
forcing of the model will cause a change in the values of a certain 
predicted variable. The difference between the perturbed and unper- 
turbed fields provides a measure of the response of the model’s atmo- 
sphere to such external forcing. In order to detect the response of the 
atmosphere caused by the different anomalies in the experiments, it is 
essential to estimate the inherent variability of the model itself. 
This is an estimate of the random variations caused by the model's 
internal, nonlinear dynamics (noise). If the signal from the perturba- 
tions or from the anomalies exceeds the noise level of the model, then 
it is considered to be a significant signal. 

For estimating the noise level of the model, a set of runs is per- 
formed to simulate the same period of integration covered by the control 
run. Each run is done by starting with different random perturbations 
in the wind and temperature fields. For each case the random initial 
perturbations have a mean of zero, and standard deviations of 1®C in the 
temperature and 4 m/s in the horizontal wind components. A statistical 
analysis is done using the results of the random initial condition runs 
to calculate the variance of the temperature and geopotential fields 
caused by the inherent variability of the model. In this study we use 
the model simulation from four different, randomly perturbed initial 
conditions, together with the control run, to give four deg^ _*es of 


freedom. 
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We define the time average of a variable at a grid point ij to 
be , and the ensemble average over five runs of this time average to 
be The geographic variance, for the variable X, can be calculated 
from the relation 


(St) 


2 

ij 



((L.) 



( 4 . 1 ) 


The standard deviation at each grid point is the positive square 
root of the above expression. The geographic distribution of the stan- 
dard deviations for the tempera ure and geopotential fields is calcu- 
lated for each of the three-month simulations. 

The difference between the monthly average anomaly case (2) and the 
monthly average control run (1) at each grid point gives the response of 
the atmosphere to the specified anomaly. This difference is given by 


M.. = |(X..)^ - (X.,),) 


( 4 . 2 ) 


The significance of the response (4^j) is calculated using the relation 


M. . 
f = -JJ 

(S,),J 


( 4 . 3 ) 


should be larger than one for a significant response of the atmo- 
sphere. That is, the absolute value of the difference caused by the 
anomalies must be larger than the inherent variability of the model at a 
certain grid point. The larger the value of tj^y l^he more confidence we 
can ii:. VC that the change is a result of the effect of the anomalies. 
The value can be large due to a large difference value or a small stan- 
dard deviation. According to Chervin and Schneider (1976), ^ was 
found sufficient for a 5% significant level, which means that the proba- 
bility of such a response arising from mere chance sampling is less 
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Table 4.1 Significance levels for a selection of critical values 
of 4 corresponding to four degrees of freedom. 


4 1 2 3 4 5 6 8 9 

significance 

level % 51.69 22.86 10.02 4.69 2.38 1.31 .77 .47 .31 


than 5%. The above table from Chervin and Schneider (1976) shovs the 
relationship between the normalized response, 4^^ j > significance 
level, using a two-sided statistical t test and assuming four degrees of 
freedom (the same number of degrees of freedom used in the p?;esent 
experiments) . 

4.4 Position and Magnitude of the Anomalies 

The experiments were designed to identify the response of the atmo- 
sphere to two main types of thermal anomalies. These are the anomalies 
in sea surface temperatures, and the anomalies in mid-tropospheric ther- 
mal heating. The anomalies are superimposed upon the original fields on 
the first day of the January (day 61 of the control run) simulation. 
However, to introduce their effect in a gradual way, as might be the 
case in the actual ocean-atmophere system, we apply a certain amplifying 
time function with an e folding time of 5 days using the relation 

(Anom)^j = (1 - e (Amplitude) , 

where V is the e-folding time. In order to isolate the effects of the 
different anomalies, the experiments were designed to test each local 
anomaly separately, as well as combined with the other anomalies, when 
appropriate. 
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4.4.’ Anomalies in Sea Surface Temperature 

Here, we primarily investigate the effect of El Nino in the eastern 
Pacific, and the cold anomalies in the Northern Pacific. However, as 
will be noted, we use super-anomalies in the Northern Pacific region. 
These are found necessary for trying to identify the atmospheric re- 
sponse in middle latitudes, where small values of SSTA and the short 
period of integration by the model (2-3 months) are not sufficient to 
detect the phase and amplitude of the response. Super-anomalies in 
those regions were used by many authors for the same reasons mentioned 
(e.g., Chervin et al . , 1980, Shukla and Bangaru, 1981). The sea surface 
temperature is maintained constant through the whole period of integra- 
tion, except in the periods of their initial gradual influence. 

Experiment 1 

This experiment describes the effect of the cold sea surface temper- 
atures in the central and western portions of the North Pacific as they 
were observed during the Fall and Winter 1976-1977. However, the actual 
values are much exaggerated in order to identify the response, for the 
reasons mentioned above. Fig. (4.1a) shows the location and magnitudes 
of the anomalies for this experiment. 

Experiment 2 

The effect of El Nino is demonstrated by this experiment. This 
warm water phenomenon lasts for 16 months on the coast of Peru, and it 
has bjen found to have significant teleconnections with middle latitude 
weather. It is important to identify its effect in midlatitudes using 
the quasi-geostrophic model and compare the model's response with the 
results of other primitive equation GCMs. Fig. (4.1b) shows SSTAs for 
Experiment 2. 
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Experiment 3 

This case represents the combined anomalies in experiments 1 and 2 
with an additional warm anomaly along the west coast ol ^ United 
States. Instead of super-anomalies in the Northern Pacific ' use more 
realistic anomalies, with a maximum amplitude of 4®C (Fig. 4.1c). 

4.4.2 Tropospheric Thermal Anomalies 

Observations in the tropical atmosphere show that the low level 
circulation over the Atlantic and Pacific is dominated by easterly flow 
on the equatorward flanks of the subtropical high pressure belts near 
30®N and 30®S. Superimposed upon this planetary-scale pattern are weak 
synoptic disturbances which move, unlike the middle latitude distur- 
bances, from east to west. These tropical disturbances are called the 
easterly waves. The cloudiness in the ITCZ is modulated by the easterly 
waves. The cloud clusters which are observed along the ITCZ are associ- 
ated with precipitation zones that accompany the regions of ascent in 
the easterly waves. The phase speed of such waves is 8-10 m sec ^ 
(Holton, 1979) and the longitudinal separation of the cloud bands is 
about 3000-4000 km, corresponding to a period of 4-5 days. 

Diagnostic studies indicate that the driving mechanism of the 
easterly waves is the release of latent heat in the convective precipi- 
tation areas accompanying the waves. Perturbations in the amplitude of 
this tropospheric thermal heating source in the tropics might cause 
certain responses in the extratropical troposphere. The monthly average 
of this enhanced heating would give an extended heat source, as often 
assumed in numerical teleconnection experiments. 

In experiments 4 and 5 we intend to investigate the effect on the 
atmospheric circulation of tropospheric thermal anomalies, which act as 
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propagating easterly waves. The anomalies will be initiated as heating 
waves at 500 mb in specific regions, using a trigonometric function in 
the form 

(Amplitude) _ cos (kx + vt), 

where (Amplitude) = 5®C/day, k is the wavenumber corresponding to 
4000 km wavelength, and V is the frequency corresponding to a 5-day 
period. 

Experiment 4 

The anomaly wave in this experiment is assumed to propagate along 
the width of the equatorial Pacific. This experiment simulates the 
tropospheric cloud pattern associated with the negative phase of the 
Southern Oscillation in which the two ITCZ's, usually located on both 
sides of the equator, merge to form a cloud cluster along the width of 
the Pacific. 

Experiment 5 

The response of the atmosphere to thermal anomalies south of the 
Himalxas are investigated. The thermal anomalies propagate at latitude 
7.7°N along the Gulf of Bengal. 
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4.1 


Sea surface temperature (°C) anomaly distribution, 
a) Experiment 1; b) Experiment 2; c) Experiment 3. 
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5. ANALYSIS OF THE RESULTS 

Since the anomalies were introduced gradually (with 5 days e- 
folding time) at the beginning of January (day 61), we decided to show 
the difference between the anomaly runs of February (day 91-120) and the 
corresponding control run. This will give the model's atmosphere suf- 
ficient time to sense and distribute the effect of the different ano- 
malies. The difference (anomaly minus control) values of the various 
fields representing the 30-day averages of zonal wind component, meridi- 
onal wind component, 750 mb temperature, 250 mb temperature, 750 mb 
geopotential, 250 mb geopotential, 500 mb pressure velocity, evaporation 
rate, and surface pressure are discussed. The significance of such dif- 
ferences is investigated using the geographic distributions of the 
standard deviations of the temperature and geopotential height fields. 
The differences in the zonal averages of the different meridional trans- 
ports and the differences of the variances resulting from stationary and 
t-ur.ol.;nt eddies also are discussed. This will clarify some mechanisms 
.iC are enhanced as a result of the various anomalies. The second 
phase of our study is to investigate the areas that retain some memory 
of the anomalies after they are gone. For this purpose, the run for the 
month of March (day 121-150) is repeated with the anomalies turned off. 
A pemory parameter is introduced to determine the magnitude of the re- 
mai response to the different anomalies. 
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Before discussing the results of the different experiments, it will 
be helpful to use linear theory to predict and compare the different 
phases of the atmospheric response to the latitudinal nositions of the 
low-level thermal anomalies. Following Hoskins and Karoly (1981), e. 
steady state linearized, iuvi&cid p plane vcrticity (Q equation may be 
written in the form 



+ p v' 


= f 


8w' 
3z • 


(5.1) 


With the came assumptions, the thermodynamic energy equation is written 
in the form 


- 30' 


30 

3y 


_L ! A 

+ V* + w 


i = (V*) 


(5.2) 


where Q represents diabatic heating. Using the thermal wind relation, 
(5.2) can be modified to yield the following equation; 


fu 


3v' 

3z 



+ w' = Q, 


(5.3) 


where N is the Brunt-Vaisala frequency. Considering the thermodynamic 
equation (5.3), three balances may arise: 


1) If the meridional advection is dominant, then 

n Q H 
£ _ j u 

^ i 3u/9z f u ’ 


(5.4) 


where is the height scale of zonal velocity, = u/(3u/3z). 
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2) If the zonal advection is dominant 
9z ” fu ’ 
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or V = 


fu ’ 


(5.5) 


where Hq is the height scale of the heat source. 

3) If the vertical advection ir dominant (this can only be true 
away from the surface), then 


w' ~ Q/N2. 


In this case the balance between the p term and the stretching 
term in the vorticity equation will require 

V' ~ fQ/ON2HQ). (5.6) 

From (5.4) and (5.5) it is clea? that the larger of the values H 

u 

and Hq will determine whether the heat source will be balanced by merid- 
ional or by zonal advection. To combine (5.4) and (5.5) in one rela- 
tion, a parameti I H is defin. 1 .n; the minimum ot and Hq. In this 
case 



(5.6) and (5.7) could be combined in a single relation by defining a new 
parameter y: 

Y = f2Q/(pN2HQH). 

This parameter differentiates between the two types of advection: 

y >> 1 low level heating will be balanced by horizontal advection 

y << 1 low level heating will be balanced by vertical advection 
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In the tropics the height scale of the heat source, .s usually 
greater than 1 km, and y is small. A heat source near urface is 
balanced by an j'pward motion (w* > 0), which implies ae creation of low 
level vorticity. Using the vorticity equation (5.1), this stietching is 
balanced by a poleward meridional motion (v* > 0). Thus, the surface 
trough must be to the *’est of the heat source. 

In middle latitudes y generally large and the heating at any 
level is balanced by horizontal advection of temperature. A heat source 
is balanced by advection of cooler air from polar regions, v* < 0, in 
the Northern Hemisphere. Thus, the trough must be to the east of the 
source. 

Similar conclusions were reached by Webster (1980) in trying to 
explain the significant resp.^nse of the midlatitude atmosphere to tropi- 
cal SSTA's, in contrast to the weak signal detected from midlatitude 
SSTA^s. The tropical SSTA initiates a direct circulation in which wars 
air rises and cold air sinks, whereas in the middle latitudes the strong 
advection will develop an indirect circulation with lisrug of cold air 
and sinking of warm air. The two reverse mechanisms are associated with 
a reverse cycle for the transformation between available potential 
energy and kinetic energy. The above explanations, although ignoring 
the nonlinear interactions, still help in understanding many of the 
results of complex models. 

The calculated average zonal wind for February is shown in Appendix 
3. The Figure will h^lp in interpreting the results of the experiments 
concerning the wind field. 
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5.1 Response of the Atmosphere tc Sea Surface Temperature Anomalies 

5.^.1 Negative Temperature Anomalies in the Northern Pacific 

Experiment 1 was designed to examine the response of the atmosphere 
to a colder than normal North Pacific Ocean. This anomaly pattern may 
have a significant effect on the weather of the United States in the 
subsequent months or seasons. The experiment was initiated using super- 
anomalies, which do not resemble lae actual observations of SSTA magni- 
tudes, However, due to the limited time of integration of the model (3 
months), we were able to identify the centers of maximum differences as 
a result of the anomalies. The results could lose much of their signi- 
ficance when using much smaller anomalies. The exaggerated anomalies 
were introduced, although with different amplitudes, in many numerical 
experiments that deal with the midlatitude SSTA’s (e.g., Chervin et al. , 
1980; Shukla and Bangaru, 1981). 

As was discussed in the previous section, a low level thermal 
cooling in midlatitudes is balanced primarily by an advection of tem- 
perature irom lower latitudes. With the meridional wind directed north- 
ward, the surface trough should be to the west of the cold anomalies. 
0i course, this simple explanation will not necessarily be true in our 
experiments, since the model design includes more complex mechanisms 
(e.g., non-linearity, mountain effects) than the simple linear theory 
assumes. 

Although the anomalies were imposed on the sea surface temperatures 
in the midlatitude North Pacific (Fig. A. la), they have a large effect 
on the 500 mb zonal winds in the polar and tropical regions. The re- 
sponse of the 500 mb zonal wind (Fig. 5.1a) shows that the major wind 
anomalies appear to the east of the surface thermal anomalies. The 
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westerlies over the west coast of the United States generally are decel- 
erated, The westerlies are enhanced in the high latitudes over central 
Canada and over the tropical Pacific (Northern Hernia here). The 
Southern Hemisphere has a negligible change as a result of the anoma- 
lies. The only substantial difference is the amplification of the 
easterlies in the tropical Pacific. The difference in the 500 mb mer- 
idional wind component (Fig. 5.1b) shows that the maximum wind anomaly 
appears to be a southward flow over the eastern United States. A maxi- 
mum northward transport is seen to the west of Canada. T^nlike the 
linear theory, a southward component is initiated directly above the 
anomaly, whereas the northward flow takes place to the east and the west 
of the anomaly. 

The difference in the 750 mb temperature fields (Fig. 5.2a) shows 
that the coldest temperatures exist directly above the anomalies, illus- 
trating the importance of the change of static stability as a result of 
cooler temperatures. The cold air mass extends over the central and 
eastern parts of the United States. A center of warm air appears over 
Greenland as a result of the northward transport over the North 
Atlantic. Although this feature of the response seems to be mainly a 
result of advection of different air masses east of the anomalies, a 
relatively strong cold air mass exists near the polar region northeast 
of Eurasia. The response in the 250 mb temperatures (Fig. 5.2b) shows 
that the difference at upper levels above the anomaly is much smaller 
than in the lower level. The positive difference over the North 
Atlantic has the same sign, but with higher magnitude, than in the lower 
level. The 750 mb temperature pattern shows that, in general, the 
anomalies tend to stabilize the atmosphere. 







5.1 The horizontal wind difference (m sec'^) at 500 mb 

between the anomaly (Experiment 1) and control run for 
February, a) zonal wind; b) meridional wind. 




5.2 


The temperature difference (°C) between the anomaly 
(Experiment 1) and control run for February. a) at 
750 mb; b) at 250 mb. 
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The 750 mb geopotential height difference (Fig, 5.3a) reveals high 
pressure over western Canada and the United States and above the area of 
the anomalies. A trough of low pressure south of Greenland dominates 
the North Atlantic. North of Eurasia a secondary high pressure region 
is found. The relatively small low pressure area over eastern Eurasia, 
west of the anomalies (in agreement with the linear theory), gives an 
indication that the Siberian high becomes weaker. This low forms mainly 
as a result of the northward transport of relatively warm air west of 
the cold anomalies. The regions of cyclonic and anticyclonic vortici- 
ties induced to the eas». of the Rockies as a result of the anomalies 
have an important influence in developing the geopotential pattern over 
the United States. The sinking motion associated with anticyclonic 
vorticity, together with the southward cold air flux, act to build and 
extend the ridge over the southern parts of the eastern United States, 
The train of waves resulting from the anomalies seems to be propagating 
towards the northeast. The quasi-barotropic character of the response 
is reflected by the 250 mb geopotential height differences (Fig. 5.3b). 
There the differences have the same sign as in the lower level, except 
over the area of the anomalies, where the 750 mb positive height anomaly 
was replaced by a low anomaly at 250 mb. 

The 500 mb pressure velocity differences (Fig. 5.4) show maxima in 
thj regions of the Rocky Mountain and Himalayas. The differences show 
an inverted omega-shaped pattern over the United States, with subsiding 
motion over the west coast and east of the continent. Over eastern 
Eurasia an anomalous rising motion is initiated. Over the North 
Atlantic the anomalous descending and ascending motions force low-level 
diver, nee, with convergence to the west and east, which enhance the 
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5.3 The geopotential height difference (m) between the 

anomaly (Fxpt'riment 1) and control run for February, 
a) at 750 mb; h) at 250 mb. 
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5.4 The vertical pressure velocity difference (10 mb 

sec at 500 mb between the anomaly (Experiment 1) and 
control run for February (negative values upward, 
positive values downward). 
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trough there over the North Atlantic. Over the westf^n North Pacific a 
descending motion is initiated above the area of the a omalicL . 

The moisture supply to the atmosphere (evaporatic rate) is largely 
decreased in the area of the cold temperatures (Fig. ".5). It is in- 
creased elsewhere in the same latitude belt over the ' cans. On a 
global basis, there is a tendency for the evaporation ra'o to decrease 
as a result of the cold SSTA. The changes are of much lower order of 
magnitude over the continental regions , 

The surface pressure differences (Fig. 5.6) show a positive area 
directly above the anomalies, which extends into the eastern parts of 
the United States. Low pressure anomalies are seen over Greenland. 
Also, noteworthy is the high pressure located in the extratropical 
latitudes of Eurasia. 

The above analysis gives a description of the main differences, 
which lasted for a sufficiently long period to appear in the monthly 
means. It is desirable to get some insight into the changes in zonal 
averages of the different transport mechanisms and in the variability of 
tne energetics resulting from both the stationary and transient eddy 
components. Table (5.1a,b) shows the difference in the zonally averaged 
transports and in the variances of the temperature and geopotential 
fields at 500 mb. 

The northward transport of sensible heat by the mean meridional 
circulation reaches its maximum in the subtropics and middle latitudes, 
compensating for the cooling effect of the anomalies. The northvard 
transport of sensible heat by the stationary eddies is dominant in he 
middle and extratropical latitudes. The transient eddies seem to be the 
major mechanism for southward sensible heat transport in the extra- 
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5.5 The surface evaporation rate difference (10 cm/day) 

between the anomaly (Experiment 1) and control run for 
February. 
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5.6 


The mean sea level pressure difference (mb) between the 
anomaly (Experiment 1) and control run for February. 



Difference (Anomaly minus Control) in the 
Zonal Average Variances and Transports 
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tropical and polar latitudes. The variance of the temperature resulting 
from stationary eddies tends to increase in all latitude belts, an 
indication of the increased land-sea temperature gradient, es;>ecially in 
the middle latitudes. This also indicates that most of the variability 
took place in the middle latitudes, as a result of the anomalies them- 
selves and of the induced circulation over land. Meanwhile, the tran- 
sient eddies seem to be damped in the latitudes north of the anomalies, 
implying less air exchange between the polar and extratropical latitude 
regions . 

The differences in the geopotential height fields show two main 
regions of mean meridional transport of potential energy: a northward 
transport in the tropical region and a southward transport in midlati- 
tudes. The transports by the stationary and transient eddies are of a 
lesser order of magnitude. Stationary eddies have a southward transport 
component north of the anomalies, where the transports by the transient 
eddies are insignificant. The differences in the variances of both the 
stationary and transient eddies are similar to those in the variances of 
the temperature fields. This result is expected since the height of a 
constant pressure level is computed using the vertical-mean temperature 
of the layer representing it. 


5.1.2 Warm Anomalies in the Equatorial Pacific 

In experiment 2, the warm SSTA’s were located in the tropical 
Pacific of both hemispheres to simulate the effect of El Nino. It is 
importar, that we carefully analyze the results of this experiment, 
con ing it with the conclusions of primitive equation general circula- 
tion models which responded to anomalies in the same location. The 
importance of such analyses stems from the use of a linear balance 


131 


0RK5INM- « 

OF POOR QUAUTT 


model, as used in the present study, which has its limitations in repre- 
senting the tropical atmosphere. The SSTA*s used here have a maximum of 
4®C and minimum of 3®C (Fig. 4.1b), closely represt ting anomalies 
observed in these regions. 

The 500 mb zonal v nd differences (Fig. 5.7a) are rr'r'^vely small. 
A positive anomaly component is seen in the tropical region of the 
Northern Hemisphere. The subtropical zonal winds a*: . generally decel- 
erated, especially over the continents. The midlatitude westerlies over 
Eurasia are enhanced. Since a considerable part of the surface anoma * 
lies lie in the Southern Hemisphere, there are some changes in zonal 
wind structure in this region, but with less magnitude than calculated 
in the Northern Hemisphere. The tropical easterJ ies and subtropical 
westerlies are enhanced. It can be noted that the anomalous westerlies 
are stronger than the easterlies in both hemispheres, which means that 
the effect of the SSTA’s is to increase the zonal momentum. This result 
agrees with maiiy primitive equation general circulation models con- 
cerning the effects of the similar anomalies (e.g., Rowntree, 19?6a). 
The meridional mass transports initiated by the anomalies (l^g. 5.7b) 
show a northward component in tn* equatorial Pacific, over the North 
American continent and over the subtropical Atlantic. The major south- 
ward mass transport takes place over the North Pacific and over the 
eastern side of Eurasia. The intersection between the northerlies and 
southerlies in the tropical Pacific lies very close to the equator, sug“ 
gesting a tendejey of the iow-ievel convergence between the northerlies 
and southerlies to shift equatorward. The order cf magnitude of the 
differences in the Nortliern Hemispliere generally is larger than that in 
the Southern Hour spiicre . 
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5.7 The horizontal wind difference (m sec ) at 500 mb 

Lween the anomaly (Experiment 2) and control run for 
oruary. a) zonal wind; b) meridional wind. 
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The differences in the temperature fields (Fig. 5.8a,b‘ reveal the 

in-situ response of the lower layer to the warm anomalies. The other 
temperature changes at remote locations are small in corparison (al- 
though they may be still statistically significant), bei ig only a frac- 
tion of a degree. The larger differences appear to be at 750 mb above 
the oceans. There are temperature differences of opposite sign in the 
Pacific to the north and south of the anomaly region. The tropical 
regions over the Atlantic are warmed, while the subtropical regions have 
a continuous area of cooler air extending from the eastern United States 
to North Africa. The relatively large differences over the oceans are a 
result of the influence of the tropical warm SSTA*s on the hydrologic 
cycle, and, therefore, on the thermal equilibrium of the atmosphere. At 
the 250 mb level the temperature differences in the South Pacific are 
almost opposite to there in the lower layer. 

The geopotential height differences at 750 mb (Fig. 5.9a) show a 
maximum in the polar region in the form of a trough of low pressure 
north of Eurasia. This feature, together with anomalously high pressure 
east of the Mediteranean Sea, develops a westerly wind component which 
intensifies the middle latitude westerlies over the Eurasian continent. 
The second feature of importance is a high pressure anomaly over the 
North Atlantic and Greenland. Most of the Pacific Ocean comes under the 
influence of a trough with two centers, the first west of Canada and the 
second in the Southern Hemispheric tropical region adjacent to the loca- 
tion of the SSTAs. The North Pacific negative height differences agree 
with results from a similar numerical experiment done jy Rowntree 
(1976a). The ' eutian low is deepened and shifted to the east, while 
a low pressure area develops near the surface anomalv location. 
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5.8 The temperature difference (®C) between the anomaly 

(Experiment 2) and control run for February. a) at 
750 mb; b) at 250 mb. 
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5.9 The geopotential height difference (m) between the 

anomaly (Experiment 2) and control run for February, 
a) at 750 mb; b) at 250 mb. 
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Our own and Rowntree^s results confirm Bjerknes’ (1969) hypothesis, 
which states that the anomalies of opposite sign over the eastern tropi- 
cal Pacific and the western Pacific appear with a negative phase of the 
Southern Oscillation. 

The 250 mb geopotential height differences (Fig. 5.9b) show the 
quasi-barotropic character of the atmosphere away from the anomalies 
(teleconnections). At the same time, directly above the anomalies, the 
negative pressure differences in the lower layer are changed to positive 
height anomalies in the upper layer, corresponding to a positive mid- 
tropospheric temperature anomaly. 

Hoskins and Karoly (1981), using a linear model, concluded that 
mid-tropospheric thermal anomalies in the subtropics tend to induce a 
wave pattern to the northeast. However, the results of the present 
experiment do not support this argument. In a nonlinear model with 
topographic forcing and differential diabatic heating, positive SSTAs in 
the equatorial Pacific are not sufficient to produce the proposed geo- 
potential pattern. 

The differences in the pressure velocity (Fig. 5.10) show strong 
rising motion in the region of the anomalies, with the maximum to the 
nortn. This also agrees with Rowntree’s experiment, where in the low- 
level convergence and ascending motion took place near the area of the 
sea surface temperature maximum. Descending motion appears east of the 
United States. Over the western Pacific, in the Northern Hemisphere, 
ascending motion takes place. 

Ai expected, the evaporation rate (Fig. 5.11) increases in the area 
of th positive SSTA. A secondary evaporation maximum appears in the 
su i^opical Atlantic in the Northern Hemisphere and in the subtropical 



-4 

5.10 The^^vertical pressure velocity difference (10 mb 
sec ) at 500 mb between the anomaly (Experiment 2) and 
control run for February (negative values upward, 
positive values downward). 
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5.11 The surface evaporation rate difference (10 cm/day) 
between the anomaly (Experiment 2) and control run for 
February. 
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Pacific in the Southern Hemisphere. The tropica^re^ift^lP^not a source 
of atmospheric moisture, due to the excess of precipitation over evapor- 
ation. The subtropical regions, which act as moisture sources, are very 
sensitive to any changes in the lower layer water cont nt and tempera- 
ture. The absolute value of the Southern Hemisphere differences gen- 
erally are greater than those in the Northern Hemisphere due to the 
increased surface temperature (summer hemisphere). Again, most of the 
differences take place over the oceans. 

The differences in the sea level pressure (Fig. 5.12) summarize the 
major features of the atmospheric response to the warn SSTA's in the 
equatorial Pacific. A surface low is induced directly above the warm 
surface anomalies. The center of the low pressure appears to the west 
of the maximum SSTA in agreement with linear theory. The Aleutian 
surface low is shifted to the east and deepened. The negative phase of 
the Southern Oscillation also is simulated as a negative surface pres- 
sure anomaly in the eastern Pacific and a positive anomaly in the Indian 
Ocean. The anomalies generally weaken the Siberian high, with the 
maximum negative values to the northeast of Eurasia. 

The differences in the zonally averaged transports and variances of 
sensible heat (table 5.2a) indicate that the stationary eddies are the 
major components in the middle and extratropical latitudes which trans- 
port sensible heat to the south. Transient eddies in the high latitudes 
have a northward transport of sensible heat, but with a 5 mailer order of 
magnitude. In the tropical region, transports by the mean meridional 
circulation have a higher order of magnitude than the eddy components. 
They are still smaller than the transports by the eddies in middle and 
higher latitudes. Although the surface anomalies are located in the 
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5.12 


The mean sea level pressure difference (mb) between the 
anomaly (Experiment 2) and control run for February. 
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Experiment 2 

Difference (Anomaly minus Control) in the 
Zonal Average Variances and Transports 



o; 


u 

CM 

3 

• 


m 

« 

w 

U 



4> 

a 

f— 4 

6 


QJ 


H 

H 



ORIGINAL PAGE IS 
OF POOR QUALITY 



uo 

os 




rH 

<N 

cn 

CM 

o 

UO 

• 

o 

Mf 

OS 

o 

00 

<r 

• 


o 

• 

'N 

00 

1 

1 

• 

1 

* 



uo 


rH 


m 

uo 

CM 

00 

VO 

CM 


o 

r-* 

vO 



VO 


• 


• 

CM 


I 

1 

• 

1 

* 


OS 

CM 




Ov 

-a- 

o 

nO* 

so 

rH 

» 

o 

00 

so 

uo 

CM 

00 

• 

• 

rH 

* 

f— ♦ 

vO 

1 

1 

• 

* 

• 


CO 





cn 

f-H 

uo 


vO 

UO 

• 

o 

VO 

rH 

VO 



* 

* 

rH 

♦ 

O 

SO 

1 

* 

• 

* 

• 







SO 


sD 

CM 

so 



O 

CM 


CO 

in 

CO 

o 

• 

O 

» 

o 

uo 

• 

1 

• 

1 

• 


vO 


UO 



ov 

CM 

CO 

o 

os 



o 

o 

o 

CM 

vO 

uo 

• 

• 

• 

• 

o 

>3* 

1 

' 

1 

* 

* 


sT 

-r 

CO 



CO 

r*i 

CM 

r-4 

c*^ 

CO 

« 

o 

O 

o 

CM 

-3* 

00 

• 


• 


o 

CO 

1 

1 

* 

1 

* 



OS 

a^ 



so 

so 

UO 

o 

UO 

rH 


CO 

O 

o 

o 

»— 4 

o 

O 


• 


O 

CO 

• 

1 

1 

1 

* 




h3- 



o 

CM 

CO ! 

o 

00 


• 


o 

o 

CM 

o 

CO 

o 1 

• 

♦ 

o 

o 

CM 

* 

1 

1 

• 

* 



r-« 




CO 

<3* 

O 

o 

so 

l-H 

* 

O 

o 

o 

O 

o 

in 

o 1 

• 

» 

o 

o 

r~( 

* 

• 

1 

■ 

• 


CM 1 



00 




uo 

t— 1 

o 

so 


O 1 

•— < 

o 

o 

o 

• 


o 

o 

• 

o 


1 

* 

* 

1 

' 


a^ 



00 



1— ( 



o 

so 

« 

i-H 

CM 


o 

r— • 

cr 

* 

o 


• 

o 


1 

* 

o 

1 

• 








c 

u 

4J 




o 

« 

c 

1 



■ *H 

C 

<L> 

c 


x: 

w 

o ^ 

•(-4 

o 


4-J 

fO 

-H u 

w u 


1 C4 


C *-H 

U 0 

C 0 

4-1 

d 

o 

«3 3 

«3 

«3 

03 CM 

«3 O 

2 

0) O 

4-^ 

S-r-i 

4J ^ 

L4 0 


s: i-i 

Wl 

HI 

w u 

H ^ 


•H 

u 

U 

0 


<U 

>»LJ ^ 

>» OJ 

>» <U 

M-4 'w' 

M-l W 

"TD 

X) U 

X w 


O 

O QJ 

3 

1—4 0 



Vi 

•H 


W (T3 

e 

E 

Vi Qj 

W *3 

• H OJ 

s- o-* 

U 

U s_^ 

QJ -H 

-a 

4-> 4J 

O Ol 

o 

o 

U TJ 

u u 

fo <u 

CU-r-i U 

CL t/3 

Cl V) 

C TJ 

d 

^ £ 

1 U3 -c QJ 

Vi OJ 

c/5 0) 

UJ 

(Q 4J 


C w 

c 

1 d -f4 


•H d 

Ui 

m u 

(T3 "O 

m 

^ >s 

L4 0> 


u Zj B 

U -U 

U T3 

CTJ L4 

(Q -H 

cu 

^ s: w 

H tjj 

H W 

> C3 

> w 


Table (5.2b) 

Geopotential 

Latitude (North) Eq. 7.7 15.3 23.0 30.6 38.3 45.9 536 61.3 68.9 76.5 84.1 
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equatorial regions » the niaximum southward sensible t transports by 
the stationary eddies are found in the extratropical latitudes. The 
same conclusion can be drawn from the differences in the variances of 
temperature resulting from stationary eddies. The extratropical lati- 
tudes reveal a maximum negative change as a result of the prevailing 
warming trend above the continental regions. The temper^' re gradient 
between land and sea is reduced. At the same time tbr ariances of 
temperature resulting from transient eddies are increased in the middle 
and extratropical latitudes, indicating an increase in the mixing of 
different air masses. 

The differences in the potential energy transport (Table 5.2b) show 
that the transport by the mean meridional circulation is dominant in all 
latitudes. The sign of transport changes between different latitude 
belts. The transports by the stationary and transient eddies appear to 
be in opposite directions. The differences of variance resulting from 
the eddies maximize in the middle and higher latitudes and resemble the 
differences in the variance of temperature. 

5.1.3 Cold and Warm Surface Temperatures in the Pacific 

Experiment 3 represents a simulation of the effects of both SSTA 
areas assumed in experiments 1 and 2 without resorting to super-cold 
anomalies in the North Pacific. However, the cold SSTA’s are still 
exaggerated to a certain extent (maximum negative anomaly 4®C)(Fig. 
4.1c). In addition, a small area of warm anomalies is located west of 
the United States. The anomalies proposed in this experiment match the 
pattern of the SSTA’s presented by Horel and Wallace (1981). According 
to their analysis, this pattern represents the first eigenvector (con- 
taining the major variance contribution, i.e. 23 percent) of Pacific sea 
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surface temperatures. It is of interest to examine the response of the 
atmosphere to these SSTA's and to see if it matches the hypothesis given 
by Horel and Wallace. The following discussion will give a good expla- 
nation of the dominant mechanisms enhancing the teleconnection between 
tropical and middle latitude regions. 

The differences in the SCO rob zonal wind (Fig. 5.13a) show that the 
major increase in the westerlies is in the Pacific sector, indicating 
that the Pacific jet is stronger and located farther to the south. 
Since the westerlies over the United States generally are decreased due 
to the superimposed negative (easterly) component. Another aiea of 
positive wind anomalies is found north of Canada, which again enhances 
an anticyclonic wind shear to the south. Both these mechanisms are 
likely to enhance high pressure conditions northwest of the United 
States and 3 low pressure trough to its southwest. Note that the major 
differences are in the western hemisphere. In the Southern Hemisphere 
the differences are of a lesser order of magnitude. The main dif- 
ferences in tile jOO mb meridiuitai wind componeta (fig. 5.13b) are lo- 
cated over the North American continent with southward component over 
eastern United States and northward component over the west. This 
feature is a result of the amplification of the Rocky Mountains of the 
anomalous easterly flow. The anomalous easterly flow turns to the south 
before it reaches the mountains. The resulting positive relative vor- 
ticity is balanced by a decrease in f (the Coriolis force) so that 
potential vorticity is conserved. As the air moves to the top of the 
mour. . 1 , the vertical column of air shrinks, and it continues to move 
equate, ard so that the decreased depth is balanced by a decrease in f. 
On the o. ■*r side of the mountains the process is reversed: The air 
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moves northward until it reaches the original latitude, where it goes 
westward* 

The differences in the 750 mb temperature pattern (Fig. 5.14a) show 
the warm and cold pools directly above the anomalies. The soMth-north 
teiaperature gradient is associated with a westerly thermal wind com- 
ponent, which enhances the jet stream in the Pacific region. Relatively 
warm air in midlatitudes extends over North America and Eurasia, whereas 
the subtropical regions there show negative temperature anomalies. The 
Southern Hemisphere reveals smaller temperature differences, except for 
the warm center in the Southern Pacific. At 250 mb (Fig. 5.14b), the 
differences are smaller than those in the lower layer, with a warm 
center above the cold anomalies and a cold center above the warm anoma- 
lies. 

The 750 mb geopotential height differences (Fig. 5.15a) reveal 
features of a blocking ridge over the United States, and low pressure 
centers on its southwestern and southeastern sides. A low pressure area 
also is found lu aortheastern rurasia. 

The 250 mb differences (Fig. 5.15b) clearly support the hypothesis 
given by Horel and Wallace (1981). We represent in Fig. 5.16 their 
Figure 11, which illustrates the hypothesized global pattern of middle- 
and upper-tropospheric geopotential height anomalies that occurs during 
an episode of warm sea surface temperatues in the equatorial Pacific 
during a Northern Hemisphere winter. The major features of our experi- 
ment and their hypothesis are a high pressure region over Canada, low 
pressure over the central Pacific, low pressure to the southeast of the 
Uni tates and high pressure over the tropical Pacific. A comparison 
with model results when each anomaly was investigated by itself 





5.15 The geopotential height difference (m) between the 
anomaly (Experiment 3) and control run for February, 
a) at 750 mb; b) at 250 mb. 
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5.16 Schematic illustration of the hypothesized global 
pattern of middle and upper tropospheric geopotential 
height anomalies (solid lines) during a Northern 
Hemisphere winter which occurs during an episode of 
warm SSTA in the equatorial Pacific (after Horel and 
Wallace, 1981). 
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(experiments 1 and 2) proves that neither is capable of producing such a 
pattern alone. Although the mid-tropospheric thermal heating in experi- 
ment 2 is larger than in the present experimen^^ (Fig. 5.8a,b and 
5.l4i*b), the response of the atmosphere in the former experiment did 
not support the results of the linear model by Hoskins and Koroly 
(1981). The negative and positive thermal anomalies to the north and 
south, are both necessary for producing a northeastly propagating geo- 
potential height anomaly pattern, which closely lesembles the obser- 
vational study by Horel and Wallace (1981). One of the effects of both 
the cold and warm anomalies acting together is the horizontal and ver- 
tical wind shear created by the atmosphere in response to the surface 
thernal anomalies. 

The differences in the 500 mb pressure velocities (Fig. 5.17) are 
large over the United States. However, a comparison with the relevant 
pattern in experiments 1 and 2 suggests that this response is a result 
of the temperature anomalies in the North Pacific, where the anomalous 
easterly component across the Rocky Mountains creates regions of cy- 
clonic and anticyclonic circulations to satisfy the potential vorticity 
conservation. 

The evaporation rate differences (Fig. 5.18) show that the warm 
anomalies cause an increase in the evaporation while the cold anomalies 
are associated with a decrease. As was the case in the previous two 
experiments, most of the variability takes place over the oceans. 

The surface pressure difference (Fig. 5.19) associated with the 
anomalies show a low pressure area in the equatorial east Pacific, which 
results from the positive SSTA underneath, and a high pressure area over 
the cold North Pacific SSTA. The subtropical high in the Northern 
Hemispheric Atlantic Ocean is intensified. 
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5.17 The vertical pressure velocity difference (10 mb 

sec~^) at 500 mb between the anomaly (Experiment 3) and 
control run for February (negative values upward, 
positive values downward). 




5.18 The surface evaporation rate difference (10 cm/day) 
between the anomaly (Experiment 3) and control run for 
February. 
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5.19 


4 


The mean sea level pressure difference (mb) between the 
anomaly (Experiment 3) and control run for February. 
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The differences in the temperature transports and variances result- 
ing from this type of SSTA distribution (table 5.3a) show that the 
sensible heat transport by the mean meridional circulation is northward 
in the middle latitudes and southward in the low and high latitudes, 
whereas the transports by the stationary and transient eddies are south- 
ward in most latitudinal zones. The southward transport by the sta- 
tionary eddies is the dominant mechanism north of 45®. The variance of 
temperature resulting from the stationary eddies is large in middle 
latitudes, unlike the results of experiment 2 where it was rather small. 
The cold anomalies in the Northern Pacific cause a larger land-sea 
temperature gradient. The variances of temperature resulting from 
transient eddies are similar to the stationary eddies; the two experi- 
ments (2 and 3) have opposite signs in middle and extratropical lati- 
tudes. Unlike experiment 2, the activity of the transient eddies is 
decreased in the present experiment. This result may give some insight 
into the variability of the atmosphere associated with the two anomaly 
patterns. Witli only a warm duomaly in the C jn itorial Pacific, the 
response of the atmosphere is transient. If the cold anomaly in the 
North Pacific is included, the response of the atmosphere is more of a 
blocking or semi-permanent type. In other words, the latter anomaly 
causes a change in the amplitude of the quasi-stationary waves and 
reduces the activity of the travelling disturbances. 

The differences in the transport of potential energy (Table 5.3b) 
by the mean meridional circulation show a maximum southward transport in 
the midlatitudes. The differences in the eddy components have a lesser 
ordc of magnitude, and they are northward in nearly all the latitude 
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belts. The differences in the variances of potential energy resulting 
from the stationary and transient eddies are similar to those found in 
the temperature field. 

5.2 Middle-Tropospheric Thermal Anomalies 

In this section we consider anomalies in the atmosphere above a 
certain location Tropospheric temperature anomalies may occur as a 
result of lower boundary temperature anomalies, as in the case of 
SSTA's, or as a result of changes in the troposphere itself. Examples 
of the latter are an anomalous release of latent heat of condensation, 
or changes iu the major atmospheric radiation absorbers like carbon 
dioxide or watei apor. The rase of troposphe'ir thermal anomaly re- 
sulting from the release of latent heat due to excessive condensation or 
precipitation is best described as a certain kind of propagating wave. 
It is unreasonable to assume that a stationary anomaly can exist in the 
troposphere for a long time. The release of latent heat will be max- 
imised in areas of ascending motion, with areas of minimum release of 
latent heat in regions of descent. Such thermal anomalies are expected 
to prevail in low latitudes (i.e., near the I7CZ) because of the large 
moisture content and convective activity in that region. 

In experiments 4 and 5 the thermal anomalies are represented by an 
easterly wave travelling along a certain latitude belt. The maximum 
amplitude of the wave is 5°C/day and the minimum is 0°C/day, resulting 
in an average of 2.5°C/day, which represents an excess precipitation 
rate of nearly 10 mm/d.ty. 
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5.2.1 Thermal Tropospheric Anoma lies Alo n g the Equatorial Pacific 

Experiment 4 represents the case of an equatorward displacement of 
the ITCZ over the Pacific in the form of a travelling easterly wave, 
with the associated anomalous release of latent heat of condensation. 

The differences in the 500 rob zonal wind field (Fig. 5.20a) reveal 
clearly that such anomalies have larger responses than the SSTA^s. 
These responses may be a result of the high amplitude used and of the 
direct and undiffused effect in the thermal energy equation. The lower 
boundary SSTA^s must be linked to the conditions of the atmosphere above 
it in order to get a large tropospheric anomaly. It also is noticed 
that the differences in the Northern Hemisphere are larger than those in 
the Southern Hemisphere. The intensification and southward shift of the 
subtropical jet stream are shown. The maximum positive anomalies take 
place over west Eurasia, increasing the westerlies there. The easter- 
lies are enhanced along the equatorial belt. Also the westerlies decel- 
erate over roost of the United States. The differences in the 500 mb 
meridional win<l elocity (F ^, .20b) show 1 u values near the 
Northern Hemispheric polar region, east of Greenland. A strong south- 
ward flow drifts towards the middle latitudes. Secondary differences 
exist in the form of a southward flux over eastern Eurasia and a north- 
ward flux over the North Pacific. The differences in the equatorial 
region are generally less than those at high latitudes. 

The differences in the 750 mb temperatures (Fig. 5.21a) show a 
strong pool of cold air east of Greenland and a weaker warm pool north- 
east of Eurasia. The major differences in the middle latitudes appear 
as warm air over most of the United States and east of Eurasia. In the 
Southern Hemisphere the response is more pronounced in the subtropics, 
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5.20 


The horizontal wind difference (m sec ) at 500 mb 
between the anomaly (Expe imerit 4) and control run for 
February, a) zonal wind; b) meridional wind. 
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with shorter wavelengths than already noted in the Northern Hemisphere. 
The 250 mb temperature differences (Fig. 5.21b) show a nearly zonal 
pattern in the Northern Hemisphere with a cold pool in the mldlatitudes. 
In the Southern Hemisphere the major positive differences are in the 
Indian Ocean with cold regions over the South African and Australian 
regions . 

The differences in the 750 mb geopotential heights (Fig. 5.22a) 
reveal that roost of the response takes place in the high latitudes of 
the Northern Hemisphere, mainly in the form of a deep trough north of 
Eurasia and a ridge over Greenland. The trough associated with the 
depression over northern Eurasia is deep ^ v reaches the western Pacific 
subtropics. The Aleutian low is enhanced over the west Pacific, while 
the Siberian high is weaker and the subtropical high is built up. 
Noteworthy are the anomalies of opposite sign over the United States. 
In the Southern Hemisphere, a low pressure area is dominant in the 
middle latitudes of the Atlantic with a high on the Antaractic region. 
The 250 mb pattern (Fig. 5.22b) shows that the major differences are 
nearly the same as in the lower layer. 

The differences in the 500 rob pressure velocities (Fig. 5.23) show 
rising motion over the equatorial Pacific region associated with the 
thermal anomalies tk^rre. In the polar region a major sinking area 
appears east of Greenland and a region of ascending motion of comparable 
magnitude o\et northern Eurasia. The mountains generate anomalous 
circulations in between these two regions. The tropical warm anomaly 
initiates a local low-level convergence accompanied by an ascending 
motion. The continuity of mass requires that air must descend in other 
regions. The descending air carries the properties of its original 
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5.2! The geopotential height difference (m) between the 

anomaly (Experiment 4) and control run for February, 
a) at 750 mb; b) at 250 mb. 
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-'i 

5.23 The_ ^vertical pressure velocity difference (10 mb 

sec ) at 500 mb between the anomaly (Experiment 4) and 
control run for February (negative values upward, 
positive values downward). 
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layer. The strong subsidence in the area of the Himalayas is very 
significant, since the subsiding air is of high momentum due to the 
existence of the subtropical jet stream there. The colder air coming 
from above is the major cause for building the ridge in the subtropics 
over Eurasia and for deepening the trough over the western Pacific. 
Descending air causes a low-level divergence, which builds the ridge, 
while the low-level convergence associated with the ascending air 
strergthens the trough to the east. In the United States the Rocky 
Mountains modify and amplify the effects of the anomalies. The mountain 
ranges produce a short-wave pattern that divides the country into two 
separate climate regions, as usually observed during anomalous winters 
in the United States. In tiio S;>athern Hemisphere the effects of the 
equatorial anomalies also are shown through the subsidence in the polar 
region. 

Fig. 5.24 shows the effects of middle- tropospheric thermaj. anoma- 
lies on the surface evaporation rate. The main differences are found in 
the subtropical oceans in both hemispheres. The circulation induced by 
the anomalies, specifically by the 750 mb temperature and horizontal 
wind anomalies, is the major reason for the evaporation anomalies, since 
the sea surface temperatures ire fixed. 

The surface pressure difference (Fig. 5.25) are similar to those of 
the geopotential heights at 750 mb. The major response occurs in the 
extratropical latitudes of the Northern Hemisphere, with high pressure 
east of Greenland, low pressure north of Eurasia and an extended trough 
over the subtropical west Pacific. The subtropical high is intensified. 
In the tropical Pacific area of the ^malous heating, a small surface 
low pressure anomaly appears. The differences usually increase with 
latitude in both hemispheres. 
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5.25 


The mean sea level pressure difference (mb) between the 
anomaly (Experiment 4) and control run for February. 
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The differences in the transports and variances (Table 5.4a) indi- 
cate that the transport of temperature by the mean meridional circula- 
tion is generaxiy southward, dominating the tropical and subtropical 
regions, whereas in the middle and higher latitudes the southward trans- 
port by the stationary eddies is a major influence. The transport of 
sensible heat by the transient eddies appears to increase poleward. The 
order of magnitude of the differences in the variances is larger than 
those calculated in the previous experiments. The variance of tempera- 
ture resulting from stationary eddies is decreased in the middle and 
high latitudes, whereas that resulting from transient eddies is in- 
creased, especially in the subtropics and the middle latitudes. 

The differences in the transport of the potential energy (Table 
5.4b) by the mean meridional circulation show a northward component at 
all latitudes except high latitudes. These estimates indicate the 
enhancement of the Hadley circulation. The differences in the variances 
of the geopotential height resembles those of the temperature fields. 

In comparison with Experiment 2 (SSTA in the equatorial Pacific), 
the pesent experiment has a similar locations of maxima and minima. The 
magnitude of the differences is larger in this experiment due to the 
larger anomaly amplitude used. 

5.2.2 Thermal Tropospheric Anomalies over the Gulf of Bengal 

In Experiment 5 the westward-propagating thermal anomaly is located 
over the Gulf of Bengal at latitude 7.7°N. The latitudinal width of the 
anomaly wave is nearly half that assumed over the equatorial Pacific 
(experiment 4). At the same time its lower boundary contains both land 
and ocean surfaces. 
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The response in the 500 mb zonal velocity (Fig, 5,26a) is shown to 
be almost zonally syimnetric. The subtropical jet stream is enhanced, 
with the maximum westerly component north of the anomaly region. The 
midlatitude westerlies are decreased in both the eastern and western 
hemispheres. The maximum differences take place in the extratropical 
latitudes, with a superimposed westerly component. In the Southern 
Hemisphere the easterlies in the tropical region and the subtropical 
westerlies are strengthened. The differences in the 500 mb meridional 
wind (Fig. 5.26b) show that the western parts of the continents have 
southward flow, whereas the eastern parts have northward flow. Again 
the maximum differences appear to be in the higher latitudes. 

The differences in the 750 mb temperatures (Fig. 5.27a) show warm 
tropical and subtropical regions in the vicinity of the anomaly. 
Northern Hemisphere middle and extratropical latitudes are characterized 
by colder air over the continents, with the major cooling over northern 
Eurasia. The Southern Hemisphere does not show any large differences 
except for the warm pool south of the Indian Ocean. The 250 mb tempera- 
ture differences (Fig. 5.27b) are mainly zonal in the two hemispheres 
with a thermal ridge in the equatorial region that extends to South 
Africa . 

The differences in the 750 mb geopotential height (Fig. 5.28a) show 
high pressure areas to the west of the continents in the Northern 
Hemisphere for high latitudes. The eastern parts of the continents 
reveal low pressure anomalies of lesser magnitude. In the lower lati- 
tudes the anomalies induce two centers of low pressure towards their 
northwest and southwest. In the Southern Hemisphere the maximum dif- 
ferences app ar as low pressure areas near Antarctica. The 250 mb 
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5.28 


The genpotential height difference (in) between the 
anomaly (Expe iuient 5) and control run for February, 
a) ac 750 ir<o; b) at 250 m'a. 
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pattern (Fig. 5 -28b) shows essentially the sane sign in the differences 
as given by the lower layer. 

The difference in the 500 nb pressure velocities (Fig. 5.29) show 
that the anonalics to the south of the Hinalayas force ascending notions 
in the whole region. Induced subsidence is created in the adjacent 
regions in the subtropics of both henispheres. The southerly flow which 
is initiated across the east-west extent of the Hinalayas acts to force 
the vertical notion. The North Atlantic is characterized by an area of 
ascending notion. 

The evaportion rate differences (Fig. 5.30) indicate that the najor 
decrease of evaporation is in the area of the Indian Ocean, south of the 
anomaly. The continents have relatively smaller differences, as was 
noted in the previous experiments. 

As in Experiment 4 the differences in the surface pressure 
(Fig. 5.31) increase towards high latitudes. The middle latitudes are 
characterized by high pressure over the continents, except over western 
Eurasia, where a trough of low pressure exists as a result of the warn 
anomalies. The Siberian high is well developed. 

The differences in the temperature transports and variances (Table 
5.5a) indicate that the mean meridional circulation dominates in the 
tropical and subtropical regions. It has northward and southward com- 
ponents to the south and north of the latitude of the thermal anomalies, 
respectively. In middle and extratropical latitudes, the stationary 
eddies have a relatively large influence, with mainly southward trans- 
port. At a certain latitudinal belt the differences in the variances of 
temperature resulting from the stationary and transient eddies have a 
similar sign. Both are negative in high latitudes. The variances 
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-4 

5.29 The vertical pressure velocity difference (10 mb 

sec at 500 mb between the anomaly (Experiment S) and 
control run for February (negative values upward, 
positive values downward). 
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“3 

The surface evaporation rate difference (10 cm/day) 
between the anonialy (Experioent 5) and control run for 
February, 
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resulting from the stationary eddies show maxima in middle latitudes. 
The transient eddy activity is increased significantly in the lati- 
tudinal zones adjacent to the anomalies. This result reflects the 
effect of the Himalayas in enhancing transient eddies caused by the 
thermal anomalies. 

The transport of potential energy (Table 5.5b) by the mern mer- 
idional circulation is northward with a maximum in the subtropical 
region north of the anomaly position. The transport is negative south 
of the anomaly. This result indicates that the Hadley circulation is 
generally enhanced as a result of the anomalies. The differences in the 
transports of potential energy by the stationary and transient eddies 
are of a lower magnitude. The differences in the variances of potential 
energy resulting from stationary waves ^>re positive and have a raaximiun 
in the high latitudes. The transient effects reveal two inver- 
sions, one in the subtropics north of the anomaly location and the other 
in the polar region. 

5.3 **Memory” of Different Anomaly Patterns in the Atmosphere 

It is difficult to estimate the way in which the atmosphere will 
"remember” certain anomalies after they are gone. The difficulty arises 
from the fact that the oceans in our model, which represent large reser- 
voirs of energy, do not respond, or change as a result of the atmo- 
spheric circulation above them. However, we can still gain an idea 
about the way in which the atmosphere recognizes the effects of the 
anomalies by itself. A memory parameter is defined to show the geo- 
graphic locations that display a sustained influence of the circulation 
caused by the anomalies. The analysis for this purpose is done in two 
phares using the March circulation (day 121-150). In the first phase 
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each experiment is conducted through March with the anomalies. In the 
second phase the experiment is repeated through the same period, but 
with a gradually decreasing anomalies. The removal of the anomalies is 
accomplished using an e* folding time of S days, which corresponds to the 
time over which the anomalies were introduced into the January circula- 
tion. 


At a certain grid point (..), we define the variable (X..). to 

ij ij i- 

represent the March anomaly run. The same variable without anosialies is 


^^ij^NA* value in the control run representing norsial March 

is (X..),. 
ij 1 

The memory parameter at grid point is given by 



(X. .)„ - (X. .)xi. 
ij'^2 ^ 1.1 ~^NA | 

(X .), - (X. .), 

xyi ^ ijn 


- 1 I X 100 . 


( 5 . 8 ) 


This parameter is calculated using the monthly average geopotential 
height of the 500 mb level at each grid point. In this way, the param- 
eter shows a high value for those locations at which geopoten^ial height 
is substantially changed by the anomalies as compared to the control 
run, while the geopotential height did not change much when the anoma- 
lies were removed. We are going to consider a significant memory to be 
above 90%, and it must cover a considerable continuous area, not just a 
single grid point. 

We will first consider the cases representing the surface tempera- 
ture anomalies (experiments 1 to 3). In experiment 1 when cold SSTA are 
located in the North Pacific, (Fig. 5.32a) the memory parameter shows 
large values over the eastern tropical Pacific. This is the area where 
the El Nino occurs. This result indicates that the atmosphere in this 
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region could preserve the effect of the SSTA to their north for longer 
time periods than could other places over either the oceans or the con- 
tinents. At the same time the differences in this region, are not large 
enough to cause a significant response of the atmosphere, as was dis- 
cussed in section 5.1.1. 

In Experiment 2, the SSTAs are located in the e.istern tropical 
Pac.fic. The response of the atmosphere to these anomalies is or a 
direct nature and has a recognizable effect on the middle and extra- 
tropical latitudes. As Figure 5.32b shows, the effect of the anomaly is 
not retained in any region after it is removed, except in the polar 
latitudes of the Southern Hemisphere. The last case in this category is 
experiment 3, where a dual-type anomaly was located in the northern and 
equatorial Pacific. The memory of the atmosphere (Fig. 5.32c) seems to 
be improved relative to the previous case. Again there are small dis- 
crete areas of high memory in the Northern Hemisphere and a large area 
in the Southern Hemisphere. An area of large memory appears west of 
Canada in high latitudes, but it is not large enough to form an air mass 
with considerable impact on in the atmosphere. 

The thenaal anomalies in the middle troposphere (experiments 4 and 
5) generally show a more continuous memory pattern (Fig. 5.32d,e), which 
spreads over large areas. We should remember the relatively large 
amplitude of the anomaly used in these experiments. In experiment 
with the thermal heat source propagating in the tropt -phere of the 
equatorial Pacific, the maximum memory regions are spre'^d over the 
subtropics north of the anomaly location, and in the middle latitudes. 
In Experiment 5 the thermal heat source is located along the Indian 
Ocean. The high memory appears to the north of the anomaly location in 


the middle latitudes of Eurasia. 
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5.4 Determination of the Noiie Level 

It is necessary to determine the noise level of the model- generated 
climatological statistics in order to judge the significance of the 
response calculated by the different experiments. The geographic dis- 
tribution of the standard deviation of temperature and geopctential 
height is calculated using the control lun and four runs with randasly 
perturbed initial conditions. The initial perturbations are iaiposed on 
the velocity and cet'perature fields at day 61 (first of January). 
Perturbed values have a zero average, with standard deviations equal to 
1°C and 4 m sec ^ for the temperature and wind, respectively (Shukla and 
Bangaru, 1980). The model runs for January and February, and their 
standard deviations are calculated using the 30 day averages of the five 
runs with their different anomaly fields. Analyses of the standard 
deviations of the temperature and geopotential height fields at 750 mb 
(Fig. 5.33a,b) show that the variability is mostly in the Northern 
Hemisphere middle and high latitudes. This is a result of the enhanced 
transient eddies during winter in these regions. The values decrease 
equatorward in both hemispheres. The maximum variability of temperature 
is between SO^N and 70°N to the west of Canada and the east of Asia. 
This feature is in qualitative agreement with a calculation dons by 
Crutcher and Meserve (1970) for the 850 mb temperatures using data from 
January of 14 years. The standard deviation of the geopotential height 
displays the same features as that of the temperature field, but shifted 
by about 90° to the east. This shift is a result of the baroclinic 
nature of the disturbances occupying these latitudes. 

As was discussed in Chapter 4, a significant atmospheric response, 
with ,'onfidence above 95%, requires the significance ratio to be > 4. A 
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5.33 Geographical distribution of the standard deviations 
(St)., from five Februaries (the normal February and 
four^^'Vandoin perturbations of the normal) simulated by 
the model . a) 750 mb temperature ('*C); b) 750 mb 
geopotential (m) . 
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coaparison between the standard deviation values and the resulting 
differences from the previous experiments shows that this condition is 
fulfilled for most of the maxisia and ainisia. This guarantees that the 
differences which result from thermal anomalies exceed (by at least four 
times) the differences caused by the random variability of the model 
itself. 


ORKalNAL PA6S 
OF POOR QUALfn- 


ORNMIALPAOer 
OP POOR QUALT 


6. SUMMARY AND CONCLUSIONS 

We intended to investigate the ataospheric response to equatorial 
and nidlatitude thermal anomalies in the form of SSTA's or diabatic 
heating anomalies in the middle troposphere. Analyses of observational 
data suggest a possible response of the atmosphere to the Pacific SSTA's 
in the equatorial and middle latitudes (e.g., Bjerknes, 1969; Namias, 
1976, 1980). Horel and Wallace 11981) showed the possible connections 
between seven warm episodes in the sea surface temperature of the equa- 
torial Pacific and the 700 mb geopotential in the Northern Hemisphere. 
Linear and nonlinear models were developed to test the possible telecon- 
nections between the atmospheric variables at different latitudes (e.g. 
Rowntree, 1972; Egger, 1977; Chervin et al., 1976; Chervin ^ al. , 1980; 
Simpson and Downey, 1975). Most of the numerical studies confirmed the 
response of the atmosphere to tropical SSTA's. Contrary to the observa- 
tional evidence, most experiments show no significant response i-f the 
atmosphere in middle latitudes, to a midlatitude SSTA of observed ampli- 
tude. In some cases super-anomalies were used (± 12®C) to obtain a 
significant atmospheric response. 

Hoskins and Karoly (1981) used a linearized steady-state five layer 
baroclinic model to study the response of the atmosphere to thenaal and 
orographic forcng. The results for the subtropical thermal forcing 
show the long waves are set up poleward and eastward of the thermal 
source, whereas shorter wavelengths are trapped equatorward of the jet. 
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The response is generally less for aidlatitude theraal forcing than for 
the case of subtropical forcing. However, the response is crucially 
dependent on the vertical distribution of the source. 

Since the atmosphere is highly nonlinear, it is difficult to pre- 
dict the correct positions of the geopotential height extremes using a 
linear model. For this purpose, a two-level, nonlinear, quasi-geostro- 
phic (linear balance) global model is developed. The dynamic structure 
of the model is basically that of Lorenz (1960). A moisture budget 
equation is used in the lower layer. In the horizontal domain a rhom- 
boidal spectral truncation is assumed, using the associated Legendre 
functions. Pressure is used as a vertical coordinate. The finite 
difference method is used for vertical approximations. 

The forcing and dissipation mechanisms adopted in the model are: 

1) Diabatic heiting due to short and longwave radiation. The 
solar energy is specified assuming an average zenith angle. 

2) Sensible heat flux between the surface and the lower layer for 
both land and sea areas. 

3) Release of latent heat of condensation in the lower layer in 
the case of near super-saturation (85%). Convective adjustment is 
used to adapt the temperature profile in order to avoid super- 
critical lapse rates. 

4) Mechanical forcing by orography as a lower boundary (1000 mb) 
condition on vertical pressure velocity. 

51 Surface friction and friction between the two layers. 

6) The sea surface temperature is specified over the oceans as- 
suming the January mean values. A steady state surface energy 
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equilibriuni is assumed over land and ice areas in order to calcu* 

late the surface temperature. 

The explicit centered scheme (Leap-frog) is used for time integra- 
tion. For computational stability the frictional and diffusion terms 
are integrated from the previous time step. A tiBie filter is used to 
couple odd and even tioie steps. The time step used is 2 hours (120 
minutes). For calculation of one month (30 days) the model requires 
about 50 seconds on the CRAY-1 machine of the National Center for 
Atmospheric Research (NCAR) . 

Integrations are performed to simulate the climate of January , 
February and March. The model starts from a hypothetical atmosphere at 
rest with a constant lapse rate. Integration is carried out for 60 days 
assuming constant insolation (first of January). From day 61 to day 
150, the January, February and March simulations are carried on with a 
daily change of insolation. This period represents the control run for 
the three months. To deduce the inherent variability of the model (the 
noise), the control run is repeated four times with different initial 
random wind and temperature perturbations. The geographic distribution 
of the standard deviations for the temperature and geopotential height 
fields are calculated using the five runs. 

Comparison of the January climate of the model with the observed 
circulation shows that the model can simulate the main teatures of the 
large scale circulation. The transient eddies are underestimated by the 
model. This limitation is due to several simplifications in designing 
the model, such as the low order truncation, the zonally averaged 
forcing of solar radiation and the time filter operator. In spite of 
these deficiencies the model is useful for climate sensitivity studies. 
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Experiments were designed to examine the response of the atmosphere 
to regional anomalies in the thermal forcing either by a lower boundary 
SSTA or by a diabatic heating rate existing in the middle troposphere. 
The main branches of the investigations are: 

1) To determine the way in which the atmosphere responds to each 
anomaly at different locations. For example, how does the atmo- 
sphere respond to a cold SSTA in the North Pacific or warm SSTA in 
the equatorial eastern Pacific, or to both anomalies existing at 
the same time? 

2) "'j estimate the response of the atmosphere to mid-tropospheric 
thermal anomalies in the tropics. 

3) To determine the atmospheric memory for the thermal anomalies 
at different locations. In other words, if the thermal anomalies 
were removed there would be different recognition of their past 
effects by the atmosphere. 

Five experiments are used to test the above arguments. Super-anomalies 
are assumed in experiment 1 (minimum -12°C) to investigate the effect of 
SSTA's in the Northern Pacific. Experiment 2 uses the usually observed 
amplitudes of SSTA to demonstrate the effect of El Nino in the equator- 
ial eastern Pacific. Experiment 3 represents the distributions of both 
of these anomalies acting together. However, the super SSTA in the 
Northern Pacific is changed to approximately observed values (minimum 
-4°C) . Experiments 4 and 5 simulate the effect of middle-tropcspheric 
anomalies resulting from the release of latent heat in a propagating 
easterly wave. In experiment 4 the anomaly wave is propagating through 
the width of the equatorial Pacific, whereas the wave is located in the 
Indian Ocean south of the Himalayas in experiment 5. 
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Experime: s with anomalies use the results of day 60 of the control 
run as an initial state. The effects of the anomalies in each region 
are io^osed gradually using a 5 day e-folding time. The insolation is 
changed daily and the model is run for January, February and March. To 
examine the way in which the model "remembers" each type of anomaly, the 
March run (days 120-150) is repeated after gradually deleting the ther- 
mal anomalies. 

The main conclusions of the model development and the experiments 
can be summarized in the following way: 

1) The nonlinear quasi-geostrophic (linear balance) model is 
useful for studying the response of the atmosphere in middle and 
high latitudes to tropical thermal anomalies. As shown in experi- 
ment 2 (SSTA in the tropical East Pacific) , the results agree with 
many primitive equation models examining the same type of anoma- 
lies . 

2) Nonlinear interactions are important mechanisms in producing 
the correct phase of the different variables. 

3) The mountains have an important influence in modifying the 
atmospheric circulation caused by the thermal anomalies. 

4) Warm SSTA’s in the equatorial Pacific induce a mid-tropo- 
spheric thermal anomaly. The corresponding response of the atmo- 
sphere does not support the linear theory argument by Hoskins and 
Karoly (1981). This anomaly pattern is not capable of producing 
the tropospheric geopotential height anomaly suggested by Horel and 
Wallace (1981). 

5) Both warm SSTA's in the tropical Pacific and cold SSTA's in 
the North Pacific, acting together (experiment 3), are essential in 
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producing the middle- and upper-tropospheric geopotential pattern 

suggested by Hoskins and Karoly (1981) and described by Horel and 

Wallace (1981), which is associated with the negative phase of the 

Southern Oscillation. The SSTA's in these regions enhance the 

thermal wind over the subtropical Pacific. With an anomalous 

easterly component over the North American continent, the northern 

and southern parts of the United States are under the influence of 

anticyclonic and cyclonic horizontal wind shears, respectively. 

6) The variances of temperature resulting from the stationary 
eddies were increased in middle latitudes in response to the SSTA's 
in the north and equatorial Pacific. The variances resulting from 
transient eddies increased in the case of equatorial SSTA’s alone. 
These results suggest that the cold SSTA’s in the North Pacific 
favor an enhancement of semi-permanent features like blocking, 
whereas the warm SSTA’s in the equatorial Pacific enhance the tran- 
sient components in middle and high latitudes. 

7) Either cold SSTA*s in the middle latitudes or warm SSTA's in 
the equatorial Pacific are capable of causing a local increase of 
the subtropical jet stream. The north-south temperature gradient 
produces a westerly thermal wind component in the lower tropo- 
sphere. 

8) The middle tropospheric thermal anomalies have a greater 
effect on the atmosphere (both local and remote) than the SSTA*s. 

9) The way in which the atmosphere "remembers" the effects of the 
anomalies after they have been removed depends on the horizontal 
and vertical locations of the anomalies. The mid-tropospheric 
thermal anomalies seem to have a more extended and organized type 
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of memory; BM>st of the "high memory" areas appear to be in the sub> 
tropical and middle latitude regions of the Northern Hemisphere. 
10) The SSTA gradient used in experiment 3 (6®C to 8“C per 40® 
latitude) seems reasonable within the capabilities of current 
climate monitoring e.g. from satellites and ships. 
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APPENDIX I 

The Solution of the Surface Thermal 
Energy Balance Equation 

The steady-state surface thermal energy balance is represented by (2.39) 

s, = S '"g - 

where 

Q = p c c , |v I (T - T ) , 
p (i* s’ g a ’ 

and Qg = L pg ^ ‘^s^'^’a^^* 

We define Ij and I 2 such that 

I, = p c c .1 V 1 
1 p d' s' 

and 

‘2 = Ps 'dl'sl- 

The above equation can be written in the form 

F(T^) = - e^BT^'' - Ij(T^ - T^) - l2(q3(Tg) - hq^O,)) 

-I(T - 271.2). (Al.l) 

g 

(Al.l) is solved for T , using Newtons iteration method. 

s 

Differentiating (Al.l) with respect to T we obtain 

O 

F' (T ) = -4e BT ^ - 
g s g 


ii - I2 - qs^V - I. 


(A1.2) 


206 




r-.’-iv ri. 

To calculate the saturation q^(T^), and its derivative, 

q' (T ), we use a formula for the saturation vapor pressure, e (Bolton, 
8 ® 

1979). This formula provides an accuracy of 0.1% in the range -30®C < 

T < 35°C. 

8 

e (T ) = 6.112 exp (17.67 (T - 273.15)/(T - 29.65)) (A1.3) 

s g 8 8 


.622 e (T ) 

= p - e^(Tg) 


(A1.4) 


Differentiating (A1.4) with respect to T , 

8 

q.pe:(TJ 




(A1.5) 


and using (A 1.3), one obtains 


,, _ 17.67 X 243.15 

s'- g-* (T - 29.65)^ • 

8 


(A1.6) 


Substituting (A1.5) into (A1.2), we arrive at 


q P X 4302.645 

F’(Tg) = - 4e^BTg3 - Ij - I2 (p-e%T ))(T -29.65)=^ ' 

S g g 


(A1.7) 


Using (Al.l) and (Al.7), the solution is convergent in the form 


V+1 V F(T ) 

T = T - S 

8 8 f'(T '')’ 

8 


(A1.8) 


where the superscripts v and v+l indicate successive i*-eration steps. 

Iteration is performed until F(T ) is less than a small, predetermined 

8 


value. 
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APPENDIX II 

Spectral Transform of (V*fVx) 


The term (V*fVx) can be expanded in the form 
V-fVx = Vf-Vx + fV^x. 

Since f = 2fip, 

+ 2Q\N^. 


(A2.1) 


If we expand x in terms of spherical harmonics defined by (2.60), then 




3Y 


n 2fi 


n ' n' 


or 


V*(fVx) = ^ 12 x™ {-(n^-DD™ y“ , -n(n+2) d"*,, 

mn n n-1 n+1 n+1^’ 

where we have used the two recurrence relations 


3Y 


m 


(1 - M^) gjS = (»tl) - „ D“^, 


and 


with 


»-in 


Y- = D"‘ , Y'" , + D™ y" , , 
^ n n+1 n+1 n n-1’ 


_m _ ,n2-m2 ^ 

*^n - 

Applying the transform operator (2.63) on (A2.3), we obtain 
(VfVx);; = t (n^-l) d" 


(A2.2) 


(A2.3) 


(A2.4) 


(A2.5) 


(A2.6) 
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APPENDIX III 

The Simulated Zonal Wind for February 



A.l 


The simulated zonal wind (m sec 


for February. 





